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ABSTRACT 


Part 1 


The aqueous solution chemistry of the methyl- 
mercury complexes of several amino acids a8 the peptide 
glutathione has been investigated by proton and carbon-13 
Magnetic resonance spectroscopy. In an amino acid such 
as valine the nature of the complex formed is pH dependent. 
At low pH protons compete with methylmercury for the 
ligand while at high pH hydroxide ions compete with the 
ligand for methylmercury. At low pH the carboxylic acid 
group is the binding site in valine while at higher pH 
the methylmercury is bound by the amino group. Formation 
constants are determined for both these complexes. Re- 
sults obtained for methionine indicate for the first time 
Significant binding of methylmercury by the thioether 
group. The methylmercury leaves the thioether group and 
Migrates to the amino group at high pH. The formation con- 
stants for the different complexes are determined. Binding 
of methylmercury by cysteine, penicillamine and glutathione 
was investigated over a Mees of solution conditions. 

In solutions containing equimolar amounts of methylmercury 
and ligand, methylmercury binds most strongly to the ion- 
ized sulfhydryl group, with no detectable dissociation of 
the complex over the pH range 0-14. Evidence is presented 
for some protonation of the methylmercury-complexed sulf- 


hydryl group of glutathione at pH < 2. In solutions con- 
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taining two moles of methylmercury: for each mole of glut- 
athione, both methylmercury cations are bound to the sulf- 
hydryl group up to pH 4, from pH 4-8, one methylmercury 
shifts from the sulfhydryl group to the amino group, while 


above pH 10 it dissociates to form-CH.HgOH. 


3 
the kinetics of the ligand exchange ,reactions 
of the methylmercury-methylamine and methylmercury- 
glutathione complexes were studied by NMR line-broadening 
techniques. Of the possible pathways by which methyl- 
mercury can exchange between the free and complexed forms 
in the methylmercury-methylamine system, first order dis- 
sociation is the only pathway contributing significantly 
to the methylmercury exchange over the pH range 2 to 5. 
The rate constant for this reaction was determined to be 
TOSS sec +, For the methylmercury-glutathione system 
the major pathway for ligand exchange depends on the 
solution pH: below pH 2 ligand exchange occurs mainly by 
means of a proton-assisted dissociation of the complex 
while above pH 2 the exchange takes place by displacement 
of complexed glutathione ee een ardor. acca glut- 
athione. The rate constants for these two reactions were 


determined to be 600 + 200 ae eee and 5.8. 1.92 108 


m™> secu: The mechanisms of the ligand exchange reactions 


are discussed. 
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Part; 2 


‘The application of carbon-13 magnetic resonance 
spectroscopy to the characterization of the acid-base 
chemistry of amino acids at the molecular level was studied 
using the amino acids cysteine and penicillamine and the 
related compound 2-aminoethanethiol. Chemical shift 
titration curves for each of the carbon atoms of these 
three compounds were obtained. From each of these curves 
macroscopic acid dissociation constants for the overlapping 
amino acid and sulfhydryl deprotonations were evaluated. 

Equations are presented which relate the ob- 
served chemical shifts to the microscopic dissociation 
constants, for such compounds, It as) shown) :that;. to eval- 
uate the microscopic dissociation constants, one needs 
to know the chemical shift of each carbon atom in each of 
the forms present. Attempts to obtain these shifts using 
chemical shift data from model compounds were unsuccessful 
presumably due to the dependence of the carbon-13 chemical 
shift on the conformation of the compound. The chemical 
shift versus pH curves of penicillamine are discussed in 


terms of the likely conformations of penicillamine. 
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PART 1 
NUCLEAR MAGNETIC RESONANCE STUDIES OF THE 


SOLUTION CHEMISTRY OF METHYLMERCURY 
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INTRODUCTION 


The chemistry of the methylmercury cation, 
CH,Hg (II), has been the subject of considerable research, 
partly due to its unique properties (1). hie isa 
metallic species with a coordination number of one (2), 
it is large, highly polarizable and the simplest soft 
acid (2,3), and its complexation reactions may proceed 
via an associative mechanism (4) in contrast to the dis- 
sociative mechanism characteristic of most other metal 
complexation reactions (5). Because of its affinity for 
sulfur, CHHg (IT) has been used as a highly selective 
reagent for protein sulfhydryl groups for some time (6). 
More recently, it has been suggested as a chemical probe 
for unpaired bases in superhelical DNA (7). 

Another dimension of interest in the chemistry 
of methylmercury has been added by the recent discovery 
that methylmercury can be formed from inorganic mercury 
by’ microbiological processes (8,9,10). Kurdangeece. al: 
(8) suggested that mercury could be alklyated by plankton 
and other marine life. In 1968 Wood and coworkers (9) 
demonstrated that both methylmercury and dimethylmercury 
are formed trom! Ho (1i)ein methanogenic bacteria. In 1971, 
these same workers elucidated reaction pathways for both 
enzymatic and non-enzymatic methylation of mercury, the 


non-enzymatic methylation proceeding via an electrophilic 
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attack by Hg(II) on methylcobalamin and suggested that any 
microorganism capable of synthesizing methylcobalamin can 
synthesize methylmercury (10). These Studies explain the 
observation that the mercury in victims of several mercury 
poisoning epidemics was in the form of methylmercury, 
even though the industrial wastes thought to be responsible 
for the mercury pollution contained inorganic mercury 
(11,12). Methylmercury has been shown to be the most 
toxic form of mercury, forming water soluble compounds 
which attack the central nervous system, essentially ir- 
reversibly. 

Since methylmercury does not exist as CHjHg” 
but rather within a coordination sphere of water as 


Gy a knowledge of its coordination chemistry is 


CH.HgOH 


3 2 
likely to be useful in understanding the behaviour of 
methylmercury in the environment and in biological systems. 
inppart...of this <thesic,. the results of nuciear 
magnetic resonance (NMR) studies of the binding of methyl- 


mercury by the nitrogen, oxygen and sulfur groups of 


selected amino acids and peptides are described. 


A. The Aqueous Solution Chemistry of Methylmercury 


The first work on the aqueous solution chemistry 
of methylmercury was reported in 1953 by Waugh, Walton 
and Laswick (14) who determined an equilibrium constant 


for the formation of methylmercuric hydroxide as described 
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by Equations (1) and (2). 


9 7 > 


CH3Hg OH. t—OH> = CHHgOH + H,0 Cr) 
[CHHgOH] 
it Pettus ese. pea (2) 


+ ER 
[CH,HgOH., ] [OH ] 


From the pH at half neutralization in the titration of 
methylmercuric hydroxide with nitric acid and perchloric 
acid, they obtained a value of 3.1 x 10° for K)- The as- 
sumption inherent in this procedure is that the only 
methylmercury species present during the titration are 

+ 


HgOH,, and CH, 


Twelve years later, Schwarzenbach and Schellenberg 


CH HgOH. 


3 
(2) showed by pH titration tnac, 11 additron co -cne 
aquated species and methylmercuric hydroxide, a dinuclear 
complex in which two methylmercury cations are attached 

to the same hydroxide ligand is also formed. The model 
developed by Schwarzenbach and Schellenberg for the aqueous 
solution chemistry of methylmercury is described by 
Equations (1) - (4). 


+ + 
HgOH, 4 CH,HgOH > (CH,Hg) 50H + HO (3) 


CH Te rs 


3 


+ 

[(CH,Hg) 50H ] 
lation S mer eneann earn gs 
-[CH,HgOH] [CH,HgOH, ] 


9 
The values they obtained for K, and K, are 2234 5x6 L0°> cand 
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The model developed by Schwarzenbach and 
Schellenberg has recently been questioned (15,16) on the 
basis of a report by Grdenic and Zado (17) that tri- 
methylmercurioxonium perchlorate, [(CH,Hg) ,0]C10,, prec- 
ipitates from concentrated neutral solutions prepared by 
titration of methylmercuric hydroxide with perchloric 
acid, and their conclusion from synthetic studies that 
methylmercuric hydroxide does not exist as a chemical com- 
pound except possibly in dilute aqueous solution as the 
dissociation product of trismethylmercurioxonium hydroxide, 
[ (CHHg) ,0] OH. These workers claimed to have prepared 
[ (CHHg) ,0] OH by a metathetical reaction between 
[(CH3Hg) ,0]C10, and KOH in methanol and concluded from 
conductance data that (CHHg) 30° is stable in alkaline 
(methanol) solution and that the compound previously 


thought to be CH,HgOH is actually [ (CH,Hg) ,0] 0H eraats 


3 
dehydration product. These results suggest that tris- 
methylmercurioxonium ion, (CHjHg) 30°, might also be present 
in aqueous solution, in which case the model of 
Schwarzenbach and Schellenberg is incomplete. In Chapter 
III, the existence of the trismethylmercurioxonium ion 


in aqueous solution is considered further and the model 


of Schwarzenbach and Schellenberg is modified to include 


this ion also: 
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B. The Coordination Chemistry of Methylmercury 


Although the parent Hg(II) ion forms both two- 
coordinate linear and four-coordinate tetrahedral comp- 
lexes (18), the predominant coordination number of mercury 
in CH Hg (IT) is two (2). It has been suggested that the 
coordination number of two is favoured for mercury in 
CH3Hg (II) due to a smaller positive charge on the mercury 
atom in complexes such as CH3HgCN as compared to Hg (CN) , 
(299% 

To determine the tendency for higher complex 
formation, Schwarzenbach and Schellenberg (2) measured 
the solubility or CH3HgI im, O<4-M KNO, and LAO hae KIS 
They observed a solubility of 1.67 x 107? moles/1l. in 


3 


0.1 M KI compared to 1.41 x 10 ~ moles/1 in 0.1 M KNO, 


and attributed the increased solubility to formation of 


3 a 
formation constant for the CHHgI, complex to be ap- 


CH,HgI From these solubility data, they calculated the 


8 ‘ 
proximately 2, as compared to 4 x 10 for the formation 
constant of CH,HgI. They also reported maximum values 


of 3 for the formation constants of CH,Hg (CN). and 


= 14 
CHHg (SR) 5 as compared to values of 1.3 x 10 for 


16 


CHLEGCN and. 1.3.x 10 for CH,HgSCH,CH,OH. These results 


3 
are consistent with a predominant coordination number of 


one for CH3Hg (IT). 


Simpson (20) observed that the formation con- 


stant for CH,HgCN did not show any systematic change when 
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the concentration of cyanide was varied over an order of 
Magnitude and also concluded higher complexes are unim- 
portant. | 

Goggin and Woodward (19) were unable to detect 


complexes of the type [cH Hg (cw) ) (871) 


, where n is 
greater than one, by Raman spectroscopy. The Raman 


spectrum of an aqueous solution of CH,HgCN containing a 


3 
three-fold excess of potassium cyanide was identical with 
that: for CH3HgCN except for an additional intense, 
strongly polarized line at 2081 em™? which is due to the 
free cyanide ion. 

Sytsma and Kline (21) considered the question 
of chelation to methylmercury. The methylmercury com- 
plexes of some potential chelating thiophenols were ex- 
amined by proton magnetic resonance spectroscopy. On 
the basis of the spin-spin coupling constants for coup- 
ling between the methyl protons and the mercury-199 
nucleus, they concluded that methylmercury is bonded only 
to the sulfur atoms of the thiophenol. 

The mercury ion found in methylmercury complexes 
is linearly coordinated by the methyl group and another 
ligand (19,22). The halide and cyanide complexes of 
methylmercury were found to possess C3, symmetry (19), 
while a mean C-Hg-S angle of 178(2)° in the D,L-penicill- 


aminato complex of methylmercury was found by X-ray 


structure determination (22). 
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The first quantitative work on the formation 
of methylmercury complexes by nitrogen-, oxygen-, and 
sulfur-containing ligands was reported by Simpson in 
1961 (20). He determined the formation constants by pH 
methods involving competition between OH and xX for 
CH,Hg (II), where X is the conjugate base of a strong 
acid, or competition between H’ and CH,Hg(II) for Bee 
where B- is the conjugate base of a weak acid, and by 
polarographic methods. The aqueous solution model he 
used was based on the model of Waugh, Walton and Laswick 
which did not include the dimeric species, (CH3Hg) OH. 
In some cases, the value obtained for the formation con- 
stant depended upon the method used. For example, the 
logarithm of the formation constant for complexation by 
the sulfhydryl group of cysteine was found to be 15.7 by 
the pH method and 16.6 by the polarographic method. The 
reason for the differences was not given. The other 
ligands studied by Simpson which are of interest in this 
thesis were acetic acid, ammonia and glutathione. He 
reported values of 3.6, 8.4, and 16.5 respectively for the 
logarithm of the methylmercury complex formation constants 
of these ligands. In his summary, Simpson states that, 
of the functional groups present in proteins, methyl- 
mercury would bind first to the sulthydry ) guoup.. ] Or ete 
other potential coordination sites, the imidazole and 


amine groups are most important; at low pH methylmercury 


‘gendined ald Gh 5 x 
bow AA | 

2 amg TERE, itd 

Sep xd etine danos 3 


ee ae his Baie 


ney 


+8: amt x08) 


atastenos. mek: tem 


\ veka ‘ones, i‘) 
i 
ay 
6 ey 
Ae “ 
i! 
ee i 
f + 
‘ 
t aie 
ae 
et hon 
a ; Cs a ogee ; an ; ; ‘paeradte r 
f yt PN, 
i A og - 


binds equally to imidazole and amine groups while at pH 
greater than 6, amine binding becomes more important 
than imidazole binding. 

In 1965, Schwarzenbach and Schellenberg reported 
the formation constants for a variety of methylmercury 
complexes (2). Some of their results are given in Table 
I. They pointed out some similarities between the methyl- 
mercury cation and the hydrogen ion. Both have a pre- 
ferred coordination number of one and form complexes 
having ligand to cation ratios greater than one only at 
very high ligand concentrations. Both form polynuclear 
complexes in which more than one cation is attached to 


the same ligand: 


. 1A r 


ve Po a2 
L rg CH3HgL = (CH3Hg) oL ae (CH,Hg) 31 


For example, when L is OH , the complexes CH,HgOH and 
(CH,Hg) ,OH” are formed and when L is S the complexes 


ob 
CH Ss, and (CH,Hg) 3S are formed. 


3HgS , (CHHg) . 
Schwarzenbach and Schellenberg point out that although 
the methylmercury cation resembles the proton in the 
simple stoichiometry of its reactions, its tendency for 
reaction with various ligands is different. Whereas the 
stability of the hydrogen halides is HF >> HCl > HBr > HI, 
the stability of the methylmercury halides is in the 


opposite order. Also the proton in general has a higher 


affinity for the oxygen atoms rather than the sulfur 
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Table I. Formation Constants for Methylmercury Complexes 


of Various Ligands 


Ligand log Key 
F_ 1.50 
Cl 5225 
Br 6.62 
Ey 8.60 
gs 2162 
CHHgS 16.3 
(CH,Hg) .S y 
HOCH.,CH.,-S™ 16.12 
HO-CH.,-CH,-S-HgCH, vs 6a27 
NH, 7.60 
NH,CH,CH.NH, 8.25 
CN 3 14.1 


2 coordinated through sulfur 


b 
coordinated through carbon 


atoms of bases, whereas methylmercury prefers coordination 
to sulfur atoms. On the basis of these relative affin- 
ities, the proton is considered a hard acid, while the 
large easily polarizable methylmercury cation is stoichi- 
ometrically the simplest soft acid (3). 

The binding of methylmercury by nucleosides 
was studied by Simpson (23). Using absorption spectro- 
scopy he determined the nitrogen sites to which methyl- 
mercury binds over the pH range 0 - 14 for the 5 naturally- 
Occurring nucleosides. More recently, Mansy, Wood, 
Sprowles, and Tobias (15) published a more complete study 
of methylmercury binding to pyrimidine nucleosides and 
nucleotides using Raman Difference Spectroscopy. They 
concluded that methylmercury binds preferentially to the 
most basic site, inthis case to one of the deprotonated 
ring nitrogens in preference to an oxygen site. 
Gruenwedel and Davidson (24) assumed methylmercury binds 
to the deprotonated nitrogen sites on the nucleotides in 
a DNA structure. Since the protons displaced by the 


CH,Hg(II) in this way are involved in Watson-Crick hydrogen 


6) 
bonding, the methylmercury causes denaturation of the 
native DNA. With the exception of Mansy, et. al. (15), 
the formation of polynuclear hydroxy complexes of methyl- 
mercury has not been taken into account in the study of 


the binding of methylmercury by nucleosides and nucleo- 


tides. Gruenwedel and Davidson assumed that at pH 5, all 
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methylmercury not complexed by nucleoside or nucleotide 
is present as methylmercuric hydroxide. The results of 
Schwarzenbach and Schellenberg indicate that this is not 
the case even at relatively low total concentrations of 
methylmercury. 

The complexation of methylmercury by several 
sulfur-containing ligands has been reported (25,22,26,28, 
29), Goggin and Woodward (25) showed that in 1.3 M aqueous 
solution prepared from methylmercuric nitrate, dimethyl- 
sulphide successfully competes with nitrate ions and 
water molecules for CHHg (II)* zoOnsito form CHHg-S(CH,),”. 
The X-ray structure of solid D,L-penicillaminatomethyl- 
mercury shows clearly that the penicillamine is ina 
zwitterion form and is bonded to the methylmercury via 
a deprotonated sulfhydryl group (22). There was no 
evidence for any intra- or intermolecular coordination 
of N or O atoms to the mercury atom. The interaction of 
methylmercury with metallothionein, a sulfhydryl-rich 
protein, has been investigated (26). Somewhat surprisingly, 
thionein has a lower affinity for methylmercury than 
cadmium, whereas the affinity of glutathione, whitch con- 
tains a single sulfhydryl group, is greater for methyl- 
mercury than for cadmium (27). 

The binding of methylmercury by the sulfhydryl 
group of N-acetylcysteine has been studied by Simpson, 


Hopkins and Hague (28) who concluded that two or more 
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methylmercury cations can bind to the sulfhydryl group 
although they did not define the solution conditions under 
which this would occur. However, Wong, et.) al. (29) 
isolated.a_ 2:1 methylmercury-penicillamine complex from 
alkaline solution and showed by X-ray crystallography 

that one methylmercury cation is bound by the sulfhydryl 
group and the other by the amino group. 

Eigen, Geier and Kruse (4) have investigated 
the dynamics of methylmercury complexation reactions by 
the temperature jump technique. They determined the 
forward and reverse rate constants for ligand displace- 


ment reactions of the type shown in Equation (5) 


HgOQH + X_ (5) 


3 


k 
CH-HoX jet OH: 4.2 CH 
2 k 


Where kinds Chu, .BY., J-4,.:CN) yy. SCNE OL SO. ac8 They con- 
cluded that the methylmercury displacement reactions 
which they studied proceed via an associative mechanism, 
that is the reacting ligand displaces the complexed ligand. 
This is different from the kinetics of most metal ion 
complexation reactions which are thought to follow a 
dissociative mechanism, that is, the rate-determining step 
is the leaving of a water molecule from the inner coor- 
dination sphere of the metal ion (5). 

It has been demonstrated by proton magnetic 
resonance (PMR) spectroscopy that the methyl group in 


methylmercury is nonlabile. The presence of satellite 
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resonances which arise from spin-spin coupling of the 
protons of methylmercury to*the 16.9% of mercury nuclei 
which has a spin of 1/2 (Hg-199) in the NMR’ spectrum 
proves that there is no dissociation of the mercury- 
carbon bond. Initially some confusion existed because 
these satellites are broadened in the iodo- and bromo- 
methylmercury complexes (30). In early NMR studies this 
broadening was interpreted in terms of exchange of the 
methyl groups; however, later work showed that the ob- 
served broadening results from quadrupole interactions 
between the iodine and bromine and the protons and not 
from methyl group exchange (31). In contrast to methyl- 
mercury, the satellites of dimethylmercury collapsed in 
the presence of ee ee trichloride due to methyl group 
exchange (32). 

In 1967, Simpson (33) measured the displacement 
of hydroxide ion from the methylmercury hydroxide complex 
by cyanide using NMR line broadening techniques. His 
results supported the conclusions concerning the associ- 
ative reaction mechanism. He also demonstrated that 
polynuclear species such as (CH3Hg) ,OH” are probably not 
of kinetic importance and predicted that the rate of 
exchange of methylmercury between unhindered sulfhydryl 
groups would be too fast to measure by NMR line broadening 
techniques. Recently Simpson, Hopkins and Hague (28) 


reported kinetic results obtained by proton magnetic 
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resonance spectroscopy for the binding of methylmercury 
chloride by N-acetyl-L-cysteine. They concluded that 
exchange ofthe peptide between the free and complexed 
forms is dominated by a pathway involving first order 
dissociation of the complex. Due to the spin-spin 
splitting of the broadened resonances, they were only 
able to estimate the rate constant for this first order 
dissociation. Their estimated value is suspect, however, 
Since the rate constant predicted for the reaction of 
methylmercury and N-acetyl-L-cysteine from their values 
and the formation constant of the methylmercury-cysteine 
complex (20) is six to seven orders of magnitude larger 
than the rate constants for a diffusion-controlled bi- 


molecular reaction. 


C. Overview of Part I 

In view of the importance of methylmercury and 
the sparsity of information about its coordination 
chemistry with amino acid and peptide ligands, it was 
decided to attempt to characterize more completely the 
binding of methylmercury by amino acids and peptides. On 
the basis of results reported for NMR studies of methyl- 
mercury compounds, nuclear magnetic resonance spectroscopy 
was chosen as the principal technique for this study. 

Both the chemical shift of the methyl protons 


and the mercury-proton spin-spin coupling constant in 
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CHHgX compounds depend on the nature of X as illustrated 
by the results Of Hatton, et. al. (30) given in Table I1. 
These results indicate that the absolute value of the 
spin-spin coupling constant decreases as the substituent 
is made more electronegative. In the case of the nitrate 
and perchlorate complexes they recognized that the species 


being observed may be CH HgOH,”. The dependence of the 


s 
spin-spin coupling constant on the nature of X is further 
illustrated by the linear eee On observed between 
spin-spin coupling constant and the aqueous solution pK, 
of the carboxylic acid in a series of methylmercury 
carboxylates in chloroform solution (34). Ina study of 
the binding of methylmercury by selected carboxylic acids 
in aqueous Coen it was found that the absolute mag- 
nitude of the spin-spin coupling constant decreases approx- 


imately linearly as both the PK, of the ligand and the 


log K, of the complex increase (35). Sytsma and Kline (21) 


F 
also observed a linear relationship between the spin- 

spin coupling constant and pK, for the ligand which they 
attributed to the fact that both the spin-spin coupling 


constant and pK, are related to the polarizability of the 


basic site to which the proton and methylmercury are bonded. 


Variation in the spin-spin coupling constant has 
also been attributed to changes in the s-character of the 
mercury-carbon bond. Federov, Kalinin, Gasanov, and 


Zakharkin (36) observed that solvation caused an increase 
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Table If. Chemical Shift of the Methyl Protons and the 
Mercury-Proton Coupling Constant of Methyl- 


Mercury in Methylmercury Complexes. 


Complex SCH, e Jou sug 
CHHgCN 0.650 178.0 
CH,HgOAc CO. 525 | 22038 
CH ,HgOH 0.442 31422 
CHHgCl 0.425 21522 
CHHgNO, O2362 227.0 
CH,HgC10, Oni7O. = 233.2 


? ppm vs. cyclohexane 
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in the spin-spin coupling constant. Goggin, Goodfellow, 
Haddock and Eary (37) also found that the spin-spin 
coupling constant depends on the ligand and increases in 
the following order: CH, < Ph << PMe, < CN = AsMe, < 
Me ONM pr. <.Cl,  < SMe, < pyridine < D,0O. They con- 


2 
cluded that the spin-spin coupling constant gives a 
measure of the metal-ligand bond strength, or at least 
the s-orbital contribution to it. The spin-spin coupling 
constant has been observed to depend linearly on the 
relative reactivities of X in Sy? type reactions and has 
been interpreted to indicate that the affinity of various 
X groups for the methylmercury cation is linearly dependent 
on the acidity of the respective acids HX (38). 

The results deeavanea above indicate that both 
the chemical shift of the methyl Deaton and the mercury- 
proton spin-spin coupling constant depend on the nature 
of the metal-ligand bond. Because of this dependence, it 
would seem that the donor group to which methylmercury is 
complexed in molecules such as amino acids and peptides 
containing more than one known potential binding site 
could possibly be identified from the magnitude of these 
parameters for the complexed methylmercury. In addition, 
the potential binding sites of the ligand can be monitored 
by means of the NMR spectra of the ligand nucle... in 
Chapter IV, the nature of the methylmercury binding by 


selected amino acids and the peptide glutathione is 
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characterized by nuclear magnetic resonance spectroscopy 
and, where possible, formation constants for the com- 
plexation have been determined. In Chapter V, the results 
of NMR studies of the dynamics of the reactions of methyl- 


mercury with several of these molecules are presented. 
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CHAPTER II 


EXPERIMENTAL 


A. Chemicals 

S-methyl-L-cysteine (Nutritional Biochemicals 
Corp.), D,L-methionine (British Drug Houses), D,L-penicill- 
amine (Aldrich Chemical Comp.), D,L-valine, tetramethyl- 
ammonium (TMA) hydroxide (Eastman Organic Chemicals), 
tertiary-butanol, methylamine, 1,4-dioxane (Fisher 
Scientific Company) were used as received. L-cysteine 
(Nutritional Biochemicals Corp.) was recrystallized from 
water as the free base. Reduced glutathione (Nutritional 
Biochemicals Corp., Terochem Laboratories, and Sigma 
Chemical Company) was washed with a water-ethanol solution 
and dried at 110°C before using. Methylmercuric hydroxide 
(Alfa Inorganics) was purified as described below. 


All other chemicals were reagent grade and were 


used without further purification. 


B. Preparation and Standardization of Methylmercuric 
Hydroxide 


The commercial methylmercuric hydroxide contained 
an acetate impurity (40), as indicated by a singlet in the 
proton magnetic resonance spectrum at 1.91 ppm downfield 
from the methyl Pee eance of sodium 2,2-dimethy1l-2-sila- 


pentane-5 sulfamic acid (DDS), and an insoluble material. 
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The methylmercuric hydroxide was purified by first passing 
an approximately 0.4 M solution in triply-distilled water 
through a 0.2 micron membrane filter two or three times 
to remove the insoluble ‘fraction.’ The filtrate was then 
passed through an anion exchange column (DOWEX 1-X8) in 
the hydroxide form to remove the acetate ions. Since 
sodium hydroxide would not interfere in the studies of 
methylmercury binding no attempt was made to remove the 
sodium hydroxide. 

The solution was then analyzed. The total 
hydroxide concentration of the solution was determined by 
adding a small excess of sodium thiosulphate to displace 
the hydroxide from the methylmercury hydroxide 


CH.,HgOH + S,0, 2 CH,HgS,0, + OH (6) 


The hydroxide was then titrated with standard nitric acid 
and the end point was located potentiometrically by 
measuring the pH as acid was added. The total hydroxide 
concentration of a typical stock solution which would be 
equal to the methylmercury concentration only if no other 
cation was present was found to be 0.366 t 0.005 M. By 
flame photometry it was found that this stock solution 
contained 0.030 * 0.005 M sodium ion while, by the NMR 
titration procedure described below, the methylmercury 
concentration of the same solution was determined to be 


0.34 + 0.01. These results indicate that, before being 
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passed through the ion exchange column, the solution con- 
sists of methylmercuric hydroxide and sodium acetate. On 
passing through the ion exchange column the sodium acetate 
is converted to sodium hydroxide. 

Ion-exchanged solutions of methylmercuric 
hydroxide were used as stock solutions and were standardized 
by one of the following two methods. The first of these 
methods is an NMR titrimetric method based on the for- 
mation of a 1:1 complex from methylmercuric hydroxide and 
thiosulphate as described by Equation (6). The titration 
was performed in an NMR tube with the titrant being added 
gravimetrically from a syringe. The endpoint in the tit- 
ration was determined by PMR by utilizing the change in 
the chemical shift of the methylmercury protons as titrant 
is added. On the PMR time scale, exchange of CH Hg (IT) 
between OH and Soca is rapid resulting in a single averaged 
resonance pattern whose chemical shift is an average of 
these parameters for the two complexes, weighted according 
to their relative concentrations. As thiosulphate titrant 
is added, the chemical shift changes linearly from that of 


CH,HgOH to that of CH,HgS,0, . The endpoint was determined 


3 
by extrapolation of the straight line portions of a plot 
of chemical shifts vs. grams thiosulfate added. From the 
grams of thiosulfate solution required to reach the endpoint 
and the densities of the solutions, the concentration of 


methylmercuric hydroxide in the stock solution was cal- 


culated. 
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Because the NMR titration procedure is time 
consuming and lacks precision, a faster and more precise 
potentiometric method based on the reaction of chloride 
ions with methylmercury to form methylmercuric chloride 
was developed. 


CHHgOH,” + Cl ¢ CH,HgCl (7) 


[CHHgC1] 
Dah SG cera ane eer oe (8) 
[CH,HgOH,”] [C17] 


The formation constant for methylmercuric chloride in 
aqueous solution is 1.78 x 10> (2). The endpoint was 
determined potentiometrically using a silver/silver 
chloride indicating electrode and a saturated calomel 
reference electrode. To minimize the competition of 
hydroxide with titrant for CH,Hg (II), the titration was 
performed in acidic solution where the major fraction of 
methylmercury exists as the hydrated cation. Initially 
the titration was attempted at pH 2 in aqueous solution, 
yielding the titration curve shown in Figure 1. The tit- 
ration curve for these conditions is drawn out, making it 
difficult to obtain precise endpoints. To increase the 
formation constant and thus the sharpness of the break at 
the endpoint, the titration was performed in a mixed water- 


ethanol solvent system. Titration curves for a range of 


solvent conditions are shown in Figure 2. In the 80% 
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Figure 1: Potentiometric titration curve for 
5.00 ml of methylmercury stock solution diluted 
to 100 ml with distilled water. The pH was ad- 
Hasted to 2 "with ‘concentrated nitric acia. * The 
titrant was 0.1024 M_ sodium chloride. 
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Figure 2: Potentiometric titration curves 
for 5.00 ml of methylmercury solution di- 
luted to 105 ml with distilled water (- e -), 
20% ethanol (- o -) or 80% ethanol (- 0 -). 
The pH was adjusted to 2 with concentrated 
nitric acid. The titrant was 0.0300 M © 


sodium chloride. 
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ethanol by volume titration the potential changed by 

0.194 v. on going from 1% before to 1% after the endpoint. 
The 80% ethanol solvent system was adopted as titration 
medium. The analysis procedure involved pipetting 5 ml 

of methylmercuric hydroxide stock solution into a 150 ml 
beaker and adding to it 15 mls of 0.3 Mimittkiel aci ditto 
reduce the pH in the final solution to approximately 2. 

80 mls of 95% ethanol were then added and the solution 
titrated with standard sodium chloride. The concentration 
of a typical stock solution was found after 4 titrations 


to be 0.4463 + 0.0005 M methylmercury. 


C. Preparation of Other Solutions 

A stock solution of tetramethylammonium nitrate, 
used as a reference for some of the chemical shift measure- 
ments, was prepared by titration of a 25% aqueous solution 
of TMA hydroxide with nitric acid to a neutral pH. 

The solutions used in the proton magnetic 
resonance (PMR) and carbon-13 magnetic resonance (CMR) 
measurements were prepared in triply-distilled water from 
pipetted amounts of stock methylmercury solution and weighed 
or pipetted amounts of the ligand. TMA nitrate, t-butanol, 
or dioxane was added for an internal chemical shift 
reference at approximate concentrations of 0.01, 0.03, and 
0.07 M respectively. To avoid dilution, the pH was ad- 


justed with concentrated nitric acid or a concentrated 
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solution of potassium hydroxide and samples of about 0.4 
ml for PMR or 2.5 ml for CMR were withdrawn at the desired 
PH values. When mole ratio measurements were being per- 
formed at constant pH, the requisite amount of ligand or 
methylmercury was added, the pH adjusted and a sample with- 
drawn. More ligand or methylmercury was added to the 
solution and the same procedure followed for all samples 
at this pH. The concentrations were corrected for the 
changing solution volume. Since high concentrations are 
necessary in the PMR and CMR experiments (greater than 
0.1 M), no attempt was made to control the ionic strength. 
When standard nitric acid was required it was 
prepared from concentrated acid and standardized with 
primary standard mercuric oxide (39). The procedure in- 
volved weighing accurately samples of dry mercuric oxide 
into Erlenmeyer flasks and then adding approximately 5 
moles of potassium iodide dissolved in a small volume of 
water for each mole of mercuric oxide. The iodide reacts 
with the mercuric oxide to form two moles of hydroxide 
for each mole of mercuric oxide according to the following 


displacement reaction. 


HgOQ + 4% + HO % HgI, + 20H (9) 


The nitric acid solution was standardized by titration of 
the hydroxide to a phenolphthalein endpoint. Standard 


solutions of sodium hydroxide were prepared from saturated 
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carbonate-free sodium hydroxide solutions and standardized 
by titration with standardized nitric acid. A stock 

solution of methylamine was prepared from a -weighed amount 
of 40% aqueous solution and standardized by titration with 


standardized nitric acid. 


D. pH Measurements 

pH measurements were made at 25 + 1°C with an 
Orion Model 801 pH meter or a Fisher Model 520 pH meter 
equipped with either a standard glass electrode and a 
fibre tip junction saturated calomel reference electrode 
or a micro-combination electrode. Saturated potassium 
acid tartrate, 0.050 M phosphate, and 0.01 M sodium tet- 
raborate solutions, DH Vaiues”’ 3.00, 4/.00) and 9.ico at 25 C 


were used to standardize the pH meter. 


E. Proton Magnetic Resonance Measurements 

Proton magnetic resonance spectra were obtained 
at 60 MHz on a Varian A-60-D spectrometer and at 100 MHz 
on a Varian HA-100 spectrometer. Both instruments were 
equipped with Varian variable temperature controllers. 
At temperatures less than 25°C, the temperature OL the 
probe was determined using the chemical shift difference 
between the methyl and hydroxyl resonances of a standard 
methanol sample and the peonenes Ae equation of the 


following set (40). 
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-98° - -48° T = 264.4 - 2.380|Aé| (10) 
-53° - -3° T = 225.4 - 2.083|Aé| (11) 
- 8° - -40° T = 195.0 = 1.810/A6| (12) 


where |A6| is the absolute value in Hz of the chemical 
shift difference. At temperatures greater than 25°C, an 
ethylene glycol sample was used with the calibration plot 
of chemical shift difference between methyl and hydroxyl 
resonances versus temperature supplied by Varian. 

Spectra were recorded at sweep rates of 0.1 
Hz/sec for the chemical shift measurements and 0.5 Hz/sec 
for the spin-spin coupling measurements. Reported data 
are the average of several scans. Chemical shifts were 
measured relative -to the central resonance of the TMA 
triplet or the t-butyl resonance of t-butanol. For all 
measurements, the chemical shifts are reported relative 
to the methyl resonance of sodium 2,2-dimethyl-2-sila- 
pentane-5-sulfonic acid (DSS) unless otherwise noted. 
Positive shifts correspond to resonances of protons less 
shielded than those of the reference. The t-butyl 
resonance of t-butanol is 1.24 ppm downfield from the 
methyl resonance of DSS while the central resonance of TMA 


is 3.17 ppm downfield from the methyl resonance of DSS. 


F. Carbon-13 Magnetic Resonance Measurements 
Car bons Nagusl. 2.420225. .— See 


Carbon-13 magnetic resonance spectra were ob- 


tained using a Bruker HFX-90 spectrometer operating at a 
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frequency of 22.63 MHz and equipped with a Nicolet 1085 
computer. The Fourier transform mode was used with proton 
decoupling. For each free induction decay signal, 8K data 
points were collected in the computer. Te eee ee were 
time averaged using 4k or 8k accumulations. The fluorine- 
19 resonance from CoFe in a coaxial capillary was used 

for the lock. The frequency range of the transformed 
Spectra was 5000 Hz. Chemical shifts are reported in 

ppm relative to the carbon-13 resonance of 1,4-dioxane 
added as an internal reference. The resonance of 1,4- 
dioxane is 67.4 ppm downfield from the TMS resonance. 
Positive chemical shifts correspond to less shielding 
than in 1,4-dioxane. The chemical shift measurements are 
eonsidered accurate to iwithin 0.05 ppm. . For all measure- 
ments the sample temperature was 25 + 1°C as measured 


with a standard thermocouple. 


G. Determination of Formation Constants by NMR Spectroscopy 
in Chapter IV ofthis thesis, Lormatien constants 
for methylmercury complexation reactions, defined by 


Equations (13) and (14), have 


+ 
CH,HgOH, + L 2% CHjHgL + H,0 (13) 


[CH.HgL] 
one (14) 


F “Gr 
[CHHgOH, ] [L] 
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been determined by NMR spectroscopy. The NMR spectrum of 
methylmercury and/or the ligand reflects the complexation, 
with the exact nature of the effect on the spectrum being 
dependent on the rate of exchange of the species being 
observed between its free and complexed forms. If the 
rate of exchange is slow, that is the reciprocal of the 
lifetime of the free and complexed forms is less than 
27A6, where Aéd represents the separation in Hz of the 
resonances of the free and complexed forms, then separate 
resonances are observed for the free and complexed forms of 
methylmercury and the ligand. The intensity of each 
resonance is proportional to the concentration of that form 
from which the formation constant can be calculated directly. 
If the rate of exchange is fast, that is the 
inverse of the lifetime is somewhat greater than 27A6, 
then exchange-averaged spectra are observed for both the 
methylmercury and ligand resonances. The observed shift 


is a weighted average of the two forms as described by 


) = P.6 + P.o (15) 


where Sop represents the observed shift, Oe represents 


S 
the shift of the complexed ligand or methylimercury, and 


6. represents the shift of the free ligand or methylmercury. 
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Equation (15), one: obtains 


6 = § 
far ORG f 
ron, ay pa (16) 


The fractional concentrations Py and Pe can be obtained 
from the observed shift and the shifts of the free and 
complexed forms, and the formation constant can be cal- 
culated from Po and Pe and the total concentration of 
ligand. 

To determine formation constants by the method 
described above, it is necessary that solution conditions 
be such that there is a significant fractional concen- 
tration of each form so that Sons differs from 5 and bu: 
This can be accomplished by utilizing competing equilibria. 
In acidic solution, competition between protons and methyl- 
mercury for the ligand can be utilized while in basic 
solutions the ligand and hydroxyl ions are competing for 
methylmercury. Using the fractional concentration deter- 
mined for the particular solution conditions from the 
observed chemical shift and the equilibrium constant for 
the competing reaction, the formation constant can then 


be calculated. 


H. Kinetic Applications of NMR Spectroscopy 
The basic theories used in the study of chemical 


exchange reactions by nuclear magnetic resonance have been 
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treated extensively in several monographs and review 
articles (41-45). The modified Bloch approach has been 
utilized in the present work to extract kinetic data from 
PMR and CMR spectra. Only a brief qualitative description 
of this approach will be given. 

Bloch (46) obtained a set of equations to des- 
cribe the time dependence of the components of the total 
nuclear magnetic moment vector per unit volume for a col- 
lection of nuclei with non-zero spins in the presence of 
a non-stationary magnetic field. The simplest chemical 
system in which chemical exchange may be studied by NMR 
is one in which the nuclei are interchanging between two 
different chemical environments, f and c, with different 
Larmor frequencies (chemical shifts). For this case, the 
time dependence of the magnetization vectors from nuclei 
in environments f£ and c can be described by Bloch equations. 
Addition of terms to these equations to account for the 
transfer of magnetization between the two sites leads to 
the modified Bloch equations (47). Solution of the coupled 
equations for the two site exchange leads to a general 


expression (Equation 17) for the intensity of the NMR 


absorption at a given frequency, ®& (in radians/sec), when 
there is exchange between the sites (48,49). in Equation 
(17), M_ represents the absorption signal, Aw, — Awe rep- 

bd 
resents the chemical shift difference of the nucleus in the 


two environments, 1/T5, and 1/T5> are the magnetic re- 
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laxation times (line widths) in the absence of chemical 
exchange, and Tie and Te represent the average time the 
nucleus spends at site c or f before an exchange occurs. 
This general equation may be simplified when the rate 
of chemical exchange satisfies certain limiting conditions. 
When the rate of chemical exchange is slow, the 
spectrum consists of two separate resonances at frequencies 
6. and Be (an. Hz~es Fors this situation the average life- 


£ 


time Tr and Ti, of the nuclei in each of the sites before 


exchange to the other site are large compared to the in- 


verse of the separation of the resonances; that is, Tp Or 


‘ -1 
T, is greater than 2m (5, 6¢) : 


determined from the extent of broadening of each of the 


These lifetimes may be 


resonances by (42) 
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' 
where Wis te 


broadened resonance of nuclei in the complexed site and 


is the width at half height of the exchange- 


gars is the width of this resonance in the. absence of 
exchange. An analogous equation may be written for the 
lifetime Tee The linewidth iy oa is related to the 
effective transverse relaxation time Too by the relation 
T = Le which takes into account the actual spin- 
spin relaxation time and broadening due to magnetic field 
inhomogeneities. 

For the condition of intermediate exchange, the 
lifetimes are of the order of the inverse of the separation 
of the resonances and the spectrum consists of partially- 
coalesced resonances. When the lifetimes are small com- 
pared to the inverse of the separation of the resonances, 
exchange is rapid and the spectrum consists of an averaged 
resonance at a frequency which is the population weighted 
average of oF and bo: [Equation GES) ]}.. “Under fast. ex= 
change conditions, the inverse of the lifetime of the com- 


plexed species, 1/T 7 is related to the various spectral 


parameters and the fractional concentrations by (50) 


2 2 
1 Pe P 4m (S¢ - 6) 
= (PW irs (Pew See 
: M1 /2708S £1 /w,b cNi/w,c 


where S¢ and bo represent the shifts of the free and 


complexed forms (in Hz), Wy /2,0BS is the observed linewidth, 


~epsifoke ath) Bo adeiad 
bas eal bakdtqmos, sal habe 


ent Od betsiat. ef at 


| ets tet ays Ra aati Syd. ailinohatade hal 
THEE, Lé auto oma sect i, ended. 
bbe st DEJ Srewesre: ot sub | 
ait l Spriahaxe- senthomebad e pres 

ire, reredae sult To amie ots to bro ote 208 


Wi/2 £ is the linewidth of the free form and W is the 
a 


172,C 
linewidth of the complexed form. 

For the methylamine system discussed in Chapter 
V exchange between free and complexed forms is rapid on 
the NMR timescale and lifetimes were calculated using 
Equation (19). For the glutathione system also discussed 
in Chapter V, the lifetimes in the individual sites were 
obtained by computing theoretical spectra as a function 


of lifetime using modified Bloch equations and then 


matching them to experimental spectra. 
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CHAPTER Cit 


METHYLMERCURY SPECIES AND EQUILIBRIA IN AQUEOUS SOLUTION 


Ast dntroduction 

The aqueous solution model used for methyl- 
mercury chemistry in the majority of studies of the 
aqueous solution chemistry of methylmercury complexes is 
that of Schwarzenbach and Schellenberg (2). This model 
was developed from pH titration data for solutions con- 
taining) 5.85) x Te) Mion ace 9Lix er M methylmercury. 
In the Schwarzenbach and Schellenberg model, CHHgOH,” 


reacts with hydroxide ion to form CH3HgOH which further 


+ : 
reacts with CH,HgOH to form (CH,Hg) .OH according to 


3 
Equations (20) - (24) 
+ ap > 
CH,HgOH, + OH 2 CH,HGOH + HO (20) 
[CH HgOH] 
Ee ree erat eer ere ee (21) 


+ -_ 
[CH,HgOH,°] [OH ] 


” “ 
CH,HgOH, + CH,HgOH 2 (CHjHg),0H” + H,0 (22) 
+ 

[(CH,Hg).OH ] 
K, = eae (23) 
H 
[CH ,HgOH,"] [CH,HgOH] 


Woodward and coworkers (51-53) identified the species 


CH,HgOH in aqueous solution by Raman spectroscopy and 


3 
assigned the band at 415 one in the Raman spectra to 
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(CHHg) ,OH”. This assignment of the Raman spectrum to 
CHHgOH has been criticized by Green (16) on the basis 

of a report by Grdenic and Zado (17) that CHHgOH does 

not exist as a chemical compound except possibly in dilute 
aqueous solution as the dissociation product of trismethyl- 
mercurioxonium hydroxide, [ (CH3Hg) 30]OH, which they 

claimed to have prepared by a metathetical reaction between 
[ (CH3Hg) 30]C10, and KOH in methanol. They concluded from 
conductance data that (CH,Hg) 40° is stable in alkaline 
(methanol) solution and thus that the compound previously 


thought to be CH,HgOH is actually [ (CH3Hg) ,0] 0H or its 


3 
dehydration product. Lorbeth and Weller (54) have shown 
that the compound reported to be [ (CH,Hg) 30] 0H by Grdenic 
and Zado is a hydrate of (CH,Hg) 40. Grdenic and Zado (17) 
also observed that a small amount of the trismethyl- 
mercurioxonium perchlorate, [(CH,Hg) ,0°]C10,", precipitated 
when a concentrated aqueous solution of methylmercuric 
hydroxide is titrated with perchloric acid, indicating 

that some trismethylmercurioxonium ion must be present in 
aqueous solution. If so an additional reaction must be 
added to the Schwarzenbach and Schellenberg model to in- 
clude formation of the species, (CH3Hg) 0°. 


+ 4 + 
CH,HgOH + (CH,Hg) 50H Jug (CH HG) Og seat Ould) 


+ 
CH,Hg),0 ] 
ae [(CH,Hg) , E RE 
[CHHgOH] [ (CH,Hg) 0H ] 
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Goggin and Woodward (53) determined by Raman 
spectroscopy that methylmercuric perchlorate is completely 
dissociated into CHHgOH,” and C10,. in aqueous solution 
as predicted by the Schwarzenbach and Schellenberg model. 
They also found that, in benzene, methylmercuric nitrate 
is not dissociated whereas in 4 M aqueous solution there 
is an equilibrium between undissociated solute and the 
3HgOH.,” and NO, . The degree of dissociation was 


£Ound to be similar to that of nitric acid -in,aqueous 


ions, CH 


solution. Clarke and Woodward (55) found by Raman spect- 
roscopy that the complex of methylmercury with methane- 
sulfonate ions is partially dissociated in 4.7 M aqueous 
solutions and completely dissociated in 1M solution. In 
the PEC Anecubeouate complex of methylmercury, the methane- 
sulfonate is characterized by monodentate bonding through 
an oxygen. The same workers concluded that methylmercury 


forms a 1:1 complex with sulfate, CH,HgSO, , which is dis- 


sociated to approximately the same extent as HSO, in 


aqueous solution. 

In this chapter, NMR and Raman results are 
reported for the aqueous solution chemistry veh e methyl- 
Mercury. These results indicate that the Schwarzenbach 
and Schellenberg model accounts for most of the methyl- 
mercury species except at high methylmercury concentrations 


where some trismethylmercurioxonium ion forms in the 


neutral pH region. This research was done in collaboration 
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with M.C. Tourangeau and Dr. C.A. Evans (56). 


B. Results and Discussion 

The PMR Bee eam of the methyl group of methyl- 
mercuric hydroxide in aqueous solution at pH 12.0 is shown 
in Figure 3. The spectrum consists of a singlet symmet- 
rically flanked by two less intense lines. The central 
resonance is due to methyl groups bonded to isotopes of 
mercury having a nuclear spin of zero while the satellites 
are due to methyl groups bonded to mercury - 199 (16% 
natural abundance, I = 1/2). The chemical shift of the 
methyl protons is given by the shift of the central 
resonance relative to the central resonance of TMA and the 
mercury-proton coupling constant is given by the separation 
between the satellite lines. The chemical shift of the 
methyl resonance and the mercury-proton coupling constant 
for a methylmercury solution containing no coordinating 
ligand other than hydroxide are pH dependent. The chemical 
shift as a function of pH is presented in Figure 4 and the 
coupling constant as a function of pH is given in Figure 
5. Both are for a solution containing 0.190 M methylmercury. 
If the model represented by Equations (20) - (25) des- 
cribes the aqueous solution chemistry of methylmercury, the 
fractional concentrations of (CHjHg) OH” and (CH Hg) ,0” 
at a given pH will be dependent on the total methylmercury 


concentration due to the unsymmetrical nature of the 
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Figure 3: PMR spectrum of the methyl group of methyl- 
mercury in a-0.190 M aqueous solution at pH = 12.0. 
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Figure 4: pH dependence of the chemical shift 
of the methyl protons of methylmercury in an 
aqueous solution containing 0.190 M methyl- 
mercury. The curves connecting the points are 
theoretical curves calculated using the con- 
stants given in the thesis. 
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Figure 5: pH dependence of the mercury-proton 
spin-spin counvling constant of methylmercury in 
an aqueous solution containing 0.190 M methyl- 
mercury. The curves connecting the points are 
the theoretical curves calculated using the con- 
stants given in the thesis. 
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equilibria for their formation. Also the species dist- 
ribution will be pH dependent, with CH 3HgOH the predom- 
inant species above a certain pH. The chemical shift of 


the protons of the methyl group and the spe 


Hq -coup Ling 
constant change. by 0.273 ppm and BynO Hz. respectively as 
the pH is increased from 0.26 to ~9 and then remain constant 
as the pH is increased further. Other measurements have 
shown that, at pH greater than 9, the chemical shift is 

the same for total methylmercury concentrations ranging 


4 


noe precipitate froma pH 12 solution’ containing 0.43°M 


Proms 0; 00437™-to 0245°M=(56)'. “Also [ (CH,Hg) ,0]C10 does 


methylmercury and 0.49 M NaClo,, indicating that the 
concentration of (CH3Hg) 0” is’ less in basic solution 
than at neutral pH. These results indicate that CH,HgOH 
is the predominant species at pH greater than 9 at these 
and lower total methylmercury concentrations. 

The non-sigmoid nature of the chemical shift 
titration curve in ‘Fig. 4 indicates that species in 


+ 
addition to CH,HgOH and CH,HgOH, are present in the pH 


3 
negion Mto.st9 TEs not possible to"observe’ the different 
species directly by NMR since time-averaged spectra are 
observed for rapidly exchanging systems such as methyl- 
mercury. However, Raman spectroscopy is capable of es- 
tablishing the presence of possible methylmercury species 


in solution if bands specific to the individual species 


are observable. 
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Raman and infrared spectra of crystalline 
[(CH3Hg) ,0]Clo, were recorded to determine the bands 
characteristic of (CH3Hg) 0° (56). Except for solvent 
bands, the Raman spectra of methanol and Wit romethane 
solutions of [(CH3Hg) ,0]Cl0, are similar to those of the 
solid, which indicates that the cation is undissociated 
in these solvents. Clarke and Woodward (57) also found 
that (CH,Hg) 0” is undissociated in acetonitrile. However, 
the spectrum of a saturated aqueous solution of 
[(CH,Hg) ,0]Cl0, shows only one intense band in the region 


uf 


500 - 600 cm, at 568 em +. In addition, some!weak to 


medium intensity bands appear in the region 400 - 500 cm +, 
These changes can be interpreted in terms of hydrolysis 
of the cation. 

Grdenic and Zado (17) had concluded from conducti- 
Metric titrations of methanolic solutions of tris (methyl- 
mercuric)oxonium salts with methanolic KOH that (CH3Hg) ,0° 
is undissociated in alkaline methanol solution. To deter- 
mine if this is the case, the chemical shift and mercury- 
proton coupling constant were measured as methanol solutions 
prepared from [(CH3Hg) ,0]C10, were titrated with methanolic 
KOH. Similar experiments were also performed in aqueous 
solution. The chemical shift results, which are presented 
in Figure 6, indicate that (CH3Hg) ,0” is not stable in 
either alkaline methanol or alkaline aqueous solution. 


The Raman results described above and the stoichiometry 
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of the titration suggest that, in aqueous solution, 


+ 
(CH3Hg) ,0 hydrolyzes to form CH,HgOH and (CHjHg) OH. 


3 
When base is added, the (CH Hg) OH” then reacts to form 
CH,HgOH. In methanol, undissociated (CH3Hg) 50° reacts 
with hydroxide, probably to form CH,HgOH and (CH3Hg) 50. 
These results are consistent with the report by Lorbeth 
and Weller (54) that the substance claimed to be 
[ (CH3Hg) ,0] 0H is actually the hydrated oxide (CH3Hg) .O°xH.,O 
The pH dependence of the species distribution 
has been determined by Raman spectroscopy by Dr. C.A. 
Evans. Raman spectra of a 0.510 M solution titrated with 
2.68 M perchloric, acid were,.recorded over the pH range 
9. -t0,.0.25, at.intervals .of. approximately 0.5 pH unit .(56). 
The high concentration was chosen to favour the formation 
of (CH3Hg) 40°. In two of the samples, at pH 6.84 and 
6.226;,,.crystals..of [ (CH3Hg) ,0]C10, slowly appeared. From 
various features in the Raman spectrum of the super- 
saturated solutions and the crystals at these pH's, the 
(CH,Hg) ,0° cation does not appear to be a major species 
in aqueous solution, even hen the solution is super- 
saturated with respect to [ (CH3Hg) ,0]C10, (56). From the 
intensity of the band at 504 cmt, which was assigned to 
methylmercuric hydroxide, K, was estimated to be 0.7 with 
a probable deviation of 0.3. Although this value of K, 
lacks precision, it:is'‘expected to be of the correct order 


of magnitude and leads to some important conclusions. 
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Calculations of species distributions, based upon this 
value for K3, assanrunction Of pH indicate: that, “ror a 
0.51 M methylmercury solution, the maximum concentration 
of (CH,Hg) ,0° is 0.019 Mand this occurs ae Srp Ores 5, 
while, at methylmercury concentrations of 0.2 and 0.05, 


tieumex tun concentr at iensuarens ex 107° 


4 


M (pH 6.4) 
and 2.0 x 10 ° M (pH 5.8). The pH dependence of the 
fractional concentrations of the various species ina 

0.2 M methylmercury solution is shown in Figure 7. 

These results are based on the model described by Equations 
(20) - (25) with log K 


= 2.31, log K. = 9.29, and 


1 2 
K, = 0.7. The Schwarzenbach and Schellenberg model 

accounts for all but a small fraction of the methylmercury 
over the pH range less than 1 to greater than 13 at the 
methylmercury concentrations which have been used in most 

of the previous studies of the solution chemistry of methyl- 


Mercury complexes (15) and those described in Chapters IV 


and V of this thesis. 
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Figure 7: Fractional concentrations of the 
methylmercury-containing species in an aqueous 
solution containing 0.200 M methylmercury. 
The fractional concentrations were calculated 
using constants given in the thesis. 
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CHAPTER IV 
The Binding of Methylmercury by Selected Amino Acids and 


by the Peptide Glutathione 


Aaa Lntreduction 

To characterize the binding of methylmercury 
by the potential methylmercury coordination sites in 
peptides and proteins, the binding of methylmercury by 
the amino acids valine, methionine, cysteine, and penicill- 
amine and the tripeptide glutathione was investigated. 
Valine was chosen as a model compound to investigate binding 
by carboxyl and amine groups; methionine, to characterize 
the interaction with the thioether group; and penicill- 
amine, cysteine and glutathione, because of the known 
affinity of methylmercury for the sulfhydryl group. The 
results of NMR studies of the complexation of methylmercury 
by these ligands in aqueous solution are presented in this 


chapter (58-60). 


B. Results 
1. The Binding of Methylmercury by Valine 
' The structural formula of valine in the 
zwitterion form is given by I. The binding of methylmercury 
by valine was investigated using the chemical shift of 


(HC) ,CHCHCO, uy 
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the methyl protons and the mercury-proton coupling con- 
stant of methylmercury. The chemical shift is presented 
as a function of pH for both a solution containing equal 
concentrations of methylmercury and valine (open points) 
and a solution of methylmercury alone (solid points) in 
Figure 8. The coupling constant data for the same solutions 
is presented in Figure 9. As shown in the figures, the 
two titration curves approach each other at both very low 
and very high pH. At very low pH, the high concentration 
of hydrogen ions results in protonation of the ligand and 
displacement of the methylmercury, while at very high pH 
hydroxyl ions compete with the ligand for the methylmercury. 
Thus complexation is favored in the intermediate pH region. 
In valine there are two potential coordination 
sites for methylmercury: the deprotonated amino and 


carboxyl groups. ,/The acid dissociation constants for 


these groups defined by Equations (26) - (29) are 
PKay = 2.29 and PKa> =) 9.810 (62) 
= + 
> 
(CH) ,CHCHCOOH + H,0 vs (CH) .CHCHCOO + H30 (26) 
+ + 
NH, NH3 
- + 
[HL] [#07] a 
K = ae 
(CH) ,CHCHCOO + HO es (CH) ,>CHCHCOO + H,0 (28) 
+ 
NH, NH, 
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Figure 8: pH dependence of the chemical shift 
of the methyl protons of methylmercury in an 
aqueous solution containing 0.190 M methyl- 
mercury (solid points) and in an aqueous solu- 
tion containing 0.265 M methylmercury and 
0.265 M valine (open points). Approximately 
one-half of the actual experimental points are 


shown. 
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Figure 9: pH dependence of the mercury-proton 
spin-spin coupling constant of methylmercury 
in anvwequeous solution containing 0, 1900M 
methylmercury (solid points) and in an aqueous 
solution containing 0.265 M methylmercury and 
0.265 M valine (open points). Approximately 
one-half of the actual experimental points are 


shown. 
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[7] [H,07] 
Ke ae aia (29) 
where PHT] represents the zwitterion of valine. 

The complexation of methylmercury by the de- 
protonated carboxyl groups of several carboxylic acids 
was studied by Libich and Rabenstein (35). These workers 
found that a carboxyl complex forms over the pH range 
1 - 9 with the amount of complex formation being strongly 
pH dependent due to protonation of the ligand at the lower 
end of this pH range and formation of methylmercuric 
hydroxide at the high pH end. It was found in a similar 
study of the complexation of methylmercury by amines that 
some methylmercury-amine complex forms over the pH range 
2 - 13 and that again the extent of complex formation is 
pH dependent due to protonation of the ligand and for- 
mation of methylmercuric hydroxide (58). The species 
distribution determined as a function of pH in the previous 
studies indicates that in molecules containing both 
carboxyl and amino groups, complexation by the carboxyl 
group will be predominant at pH less than 3 while at pH 
greater than 8, binding to the amino group will be most 
important. The NMR titration curves presented in Figures 
8 and 9 above pH 8 are very similar to the NMR titration 
curves for the binding of methylmercury by amines which 
indicates that, above this pH, methylmercury is binding 


to the amino group. 
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The binding of methylmercury by the amino group 
is described by Equations (30) and (31). The formation 


constant for this reaction 


Bk ae = 
(CH. ) ,CHCHCOO Br CH.Hg < (CH) ,CHCHCOO (30) 
NH, sik 
Hg 
3 , CH. 
H,HgL 
N 


[CH,Hg™] [L7] 


was determined in the general manner outlined in Chapter 
II Section G using chemical shift data of the type 
presented in Figure 8. To ensure that only the amino 
complex is present, data at pH greater than 9 was used. 
The fractional concentration of amine-complexed methyl- 


Mercury is given by Equation (32). 


pa, Saonsiiees (32) 
re. os 
c £ 


where Po is the fractional concentration of complexed 
methylmercury, o- is the chemical shift of the complexed 
methylmercury, Fe is the shift of free methylmercury, and 
Sops is the observed chemical shift. In this pH range, 


the "free" methylmercury is essentially all methylmercuric 


hydroxide whose chemical shift is equal to 2.35 ppm. The 
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chemical shift of the complexed methylmercury cannot be 
determined directly because at no pH is all the methyl- 
mercury in the complexed form. Consequently an iterative 
procedure was used to simultaneously evaluate ore and the 
LormMeacton constan.,. initially a chemical ‘shift,o£ : 2.27 
ppm was used, which corresponds to the shift of the nearly 
level portion of the NMR titration curve between pH 7 and 
eh ahs P. was calculated from Equation (32) using this ini- 
tial value for oe for’ the: data at pH greater than 10. Py 
was then calculated for each data point from the Equation 


POS i= Poe The concentration of the complex at each pH 


£ 


was then calculated from the relation 


[CH,HgL] = PM, (33). 


where Mi is the total methylmercury concentration. At 
PH above 10, essentially all of the "free" methylmercury 
is in the form of the hydroxide and [CH,HgOH] = PM. 
The concentration of methylmercury in the cationic form 


is calculated using Ky the formation constant of the 


hydroxyl complex. (See Section B of Chapter III). 


i [CHHgOH] 
ere Ogee er 0 ae (34) 
(2.34 x 107) 1OHM} 


Since the total methylmercury and valine concentrations 
are equal, the fractional concentration of free valine 


is also l - Pas The concentration of amino-deprotonated 
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valine is given by 


[L ] = ao,P.M (35) 


+ ! 
where a, equals Kyo/Ky + [H ] (62). In this manner a 


value was obtained for the formation constant KN from 


each data point. Using the average of the individual 
formation constants and the - used in their calculation 
a chemical shift titration curve was calculated. On was 
then changed slightly and a new set of formation constants 


obtained. The - and log Ky which resulted in the smallest 


standard deviation between calculated and observed chemical 
shift titration curves are 2.26 ppm and 7.41. The stan- 


dard deviation of log Ke is 0.08: %¥The formation constants. 


calculated from individual data points using the value 


2.26 ppm for the chemical shift, Our are given in Table III. 


The formation constant for the complexation of 
methylmercury by the carboxylic acid group of zwitterionic 
valine as defined by Equations (36) and (39) was evaluated 


from chemical shift data 


Hg’ 2%  (CH,) ,CHCHCOOHgCH (36) 


(CH) ,CHCHCOO VCH 3 


3 


+ + 
NH, NH 


[HLHgCH ,” ] 
Ree aie eee ee (37) 
[HL ] [CH,Hg ] 


57 


= 


(ee). 


6 Leen eins a 
mont o passage 

=) being fret pid ¥, 
ngirelvotee seats nk & 
Pua a _hveainaite 


- 


neat tae 
od? 


foal tina per as 1 heal 
tneabete ‘gaia alin nae & ‘ 

“ove ott EDS bis : 
sennseno> soissaaat ty , 
ee olsen 2 tay 
%0 node Same 


bedeut ave por a ts: ren 


be ab aii nae 


gs 
yin ; uy 
= =" be ' ‘ Ai, Re K 
! = i 1 fice i”, 
rt i “ ‘ 
rae oe 
S } eS 
a nce 
7 Os 
\ if i} 
' o 
T ge } i 
my 
’ ah 
' 
(Ve) . | 
; , 
, ; r 
a i ame Pe ' 
= § ie : w ‘ ' 
F “tl 
" =f { : 2 oy, ‘a ' Paes, Fae wy i) om 
ya ‘ } by ; iy Barry i : Le et yf 
i ( joy 


Table III. Calculated Values for Ke 


pH log Ky 
10.50 TAL 
10.80 7.41 
11.00 7.40 
11.30 7.43 
11.50 7.40 
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anepn Tess thans3..) The calculation of ge is complicated 


by the fact that, in (this pH region, the "free" methyl- 
mercury is distributed among several forms and, ata 
given pH, the species distribution depends on the "free" 
methylmercury concentration. Consequently, 5. is the 


weighted average of the chemical shifts of the various 


forms of "free" methylmercury. Also, as in the determination 


of ray by cannot be measured directly. he was evaluated 
by an iterative procedure similar to that used in the 


evaluation of me with the addition of an iterative 


Tee 
procedure for the evaluation of b- for the specific "free" 
methylmercury concentration. In the iterative procedure 


for evaluating 6 the concentrations of the various 


£! 
"free" methylmercury species were calculated for the "free" 
methylmercury concentration predicted initially using the 
value of O- at cheigiven pl stom tire }0e90 M methy mercury 
SOlUtIOn- (Pagures4) SnvEquation W(o2)... “From these concen- 
trations a new value of Ob is obtained which is then used 
to predict a new value for the "free" methylmercury con- 
centration. The iterations were continued until the "free" 
methylmercury concentration remained constant. In each 
iteration, the concentration Of the cationic! form of 


methylmercury for a particular "free" methylmercury con- 


centration was calculated by 


) 


2Nue ie ; (1-4 K 2 
[MeHg ] - (MeHg) 


1. 


+ free. : (38) 


[H ] fH) 
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which is derived from Equations (20) and (22). The con- 
centrations of (CH3Hg) OH” and CH ,HgOH were then calculated 
using Equations (21) and (23) and the concentration of 

the cationic species obtained from Equation (38). 
Typically, after 4 iterations, successive Py values dif- 
fered by less than 0.1 %. The concentration of carboxyl 
deprotonated valine used in the evaluation of pa was 


F 


calculated from the relation 


— 
(Hi ~] = a5P, T (39) 


where a a3 [H"] (62) and Mn is the total methyl- 


ieee Was 
mercury concentration. The values of oe and log Ke which 
resulted in the smallest standard deviation between cal- 
culated and observed chemical shift titration curves were 
2.16 ppm and. 2:'7:. The standard deviation of log Ke is 
0.2. The individual results obtained are given in Table 
Iv. The relative lack of precision in the carboxyl 
formation constant determination, as compared to the amino 
formation constant determination, arises from the fact 
that the observed shift of the methylmercury does not 
differ greatly from the shift of the free methylmercury 

in the pH region 1-5 whereas in the high pH region, where 


the amino complex forms, the difference between the 


observed and free shifts is quite large. 
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Calculated Values of K? 
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2. The Binding of Methylmercury by Methionine 
The structural formula of methionine in the 
zwitterionic form is 


H,CSCH,CH.CHCO, EL 


+ 


NH, 


Natusch and Porter found by PMR that at low pH, Hg(II) 
was bound by the thioether group (63,64) which suggests 
that the thioether group might also bind methylmercury. 
The -CH,CH.CH- protons of methionine form an AMIN, spin 
system and the methine proton gives rise to an approximate 
triplet. The shift of the central peak of this triplet 
was used to monitor interactions at the carboxylic acid 
and amino groups. While this may not represent the exact 
shift of the methine proton, it is a parameter which is 
sensitive to changes in the electronic environment due to 
deprotonation or complexation reactions. The methyl protons 
give rise-to a sharp singlet. 

The curves through the solid points in Figure 10 
represent the shifts of the methine and methyl protons 
inva 03450 M methionine solution, while the curves through 
the open points represent the shifts of these protons for 
a solution containing 0.160 M methylmercury. The curves 
through the solid points in Figure ll and 12 represent 


the shift and coupling constant of methylmercury in a 


solution containing 0.190 M methylmercury while the curves 
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CHEMICAL SHIFT, ppm vs DSS 
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Figure 10: pH dependence of the methyl 
(upper curves) and methine (lower curves) 
protons of methionine in aqueous solu- 
tions containing 0.150 M methionine 
(closed points) and 0.160 M methionine 


and 0.160 M methylmercury (open points). 
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CHEMICAL SHIFT, ppm vs TMA 


2.38 


2.30 


2.22 


2.14 


206 


0 4 8 12 
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Figure 11: pH dependence of the chemical 
shift of the methyl protons of methyl- 
mercury in an aqueous solution containing 
0.190 M methylmercury (solid points) and 
in an aqueous solution containing 0.160 M 
methylmercury and 0.160 M methionine ~— 
(open points). i 
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COUPLING CONSTANT, Hz 


200 


Figure 12: pH dependence of the methyl- 
mercury mercury proton Spin=spin Coupling 
constant in aqueous solutions containing 
0.190 M methylmercury (solid points) “and 


O60 Mime thy Imercury and +0). 160 °M métiio= 
nine (open points). 


65 


tig a i at map 
te =e 


- 


a 


ar Na pg PEG Geman 
=e = = i ani make 


through the open points represent the shift and coupling 
constant data for the solution containing 0.160 M methyl- 
mercury and 0.160 M methionine. 

It is well known that the shifts of carbon- 
bonded protons close to acidic functional groups are sen- 


Sitive to the deprotonation of that functional group (65). 


Since the methine proton is close to both functional groups, 


its chemical shift, as shown in Figure 10 for the 0.150 M 
methionine solution, moves upfield as the carboxyl group 
is deprotonated in the pH range 0-4, remains relatively 
constant in the pH range 4-7, and then shifts further up- 
field in the pH region 7-11 as the amino group is de- 
protonated. In the pH region 0-4 the observed shift is 
the weighted average of the shifts of the forms in which 


the carboxyl is protonated and deprotonated as described 


by 
6 PG rt Re (40) 
where the subscripts p and d refer to the protonated and 


deprotonated forms. Substitution of Ls + Pa = 1 into 


Equation (40) results in the following Equations for 


and Pac 
6 - 6 
Soe OBS = d (41) 
p os 54 
6 - 6 
OBS p (42) 
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Kay was then calculated from each of the observed shifts 


in this pH range using Equation (43). Similarly the 


ialBe 
sy. WO aaa Eine (43) 
p 


observed chemical shifts in the pH range 7-11 are the 
weighted average of the shifts of amino protonated and 
deprotonated forms from which the acid dissociation con- 
stant K for the amino group was calculated. The values 


A2 


for pK and pK, are 2.0 and 9.2, both with a standard 


Al 2 
deviation of 0.1. The shift of the methyl protons in 
methionine is essentially independent of pH over the pH 
range 0-14 which indicates that it is not affected by 

the state of protonation of the distant amino and carboxyl 
groups and that the thioether group is not protonated even 
under very acidic conditions. 

At pH less than 2 the chemical shift of the 
methine proton is not affected by the presence of methyl- 
mercury (see Figure 10), whereas the methyl protons are 
shifted downfield which indiicateemenae, in this pH range, 
methylmercury is bonded to the thioether group. As the 
pH is increased from pH 2 the chemical shift curves given 
in Figure 10 for the methine proton in the presence and 
absence of methylmercury become increasingly different, 
while the curves for the methyl protons approach each 


other indicating that the binding site of the methyl- 
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mercury is shifting away from the thioether group to the 
other end of the molecule. At pH ~8-9, coordination is 
exclusively to the amino end. As the pH is increased 
above 9, the complex begins to dissociate as indicated by 
the shift of the methine resonance of the methylmercury- 
containing solution towards that of free methionine. The 
overlap of the two curves at pH greater than 13.5 in- 
dicates complete dissociation of the complex at this pH. 


The chemical shift and the 1,199 


Hg coupling 
constant of the methyl protons of methylmercury presented 
in Figures 11 and 12 also support these conclusions con- 
cerning complexation of methylmercury by methionine. Thus 
the NMR data indicates that a thioether complex forms 
according to Equation (44) in the low pH region. 
CH,SCH.CH.CHCO.H + CH Hg” < CH,SCH,CH,CHCO 5H (44) 


+ ‘, 
NH, +HgCH, NH, 


The formation constant for this complex is 


2+ 
alee oe 
[HL"] [CH,Hg"] 


[H.,LHgCH 
(45) 


A 
I 


At higher pH, the methylmercury shifts to the other end 
of the molecule and potentially could be coordinated to 
either the deprotonated amino or the carboxylic acid 
groups. By analogy with valine at pH greater than: 9, the 


methylmercury is assumed to be bound by the amino group 
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as described by reaction (46). The formation constant 


S a = 
CH,SCH,CH,CHCO., + CHHg we CH,SCH.CH.CHCO., (46) 
iu hi 
H. NH»-HgCH, 
+ 
for this reaction is given by Equation (47). At inter- 


mediate pH values a complex may form in which the methyl- 


[CH ,HgL"] 
Keh@ sched) “ste =: (47) 
[L] [CH Hg") 


mercury is bonded to the deprotonated carboxyl group, while 
the amino group is protonated, although it was not possible 
to identify such a complex in the present work. 

The formation constants for the reactions des- 
cribed by actos Gas (44) and (46) were determined from NMR 
data. The formation constant for the binding of methyl- 
mercury by the thioether group was determined from the 
chemical shift of the methyl protons of methionine. The 
chemical shift was measured as a function of the ratio of 
methionine to methylmercury at pH 0.5 and pH 1.5. For this 
study each sample was prepared separately with the exact 
stoichiometry achieved by pipeting methylmercury stock 
solution and adding a weighed amount of methionine. Only 
the results at pH 0.5 were ultimately used since preliminary 
calculations indicated the possibility of significant 


carboxyl complexation at pH 1.5. The data at pH 0.5 are 
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presented in Table V. At this low pH, the amino and 
carboxylate dentates are protonated, as indicated by the 
shift of the methine proton, and "free" methylmercury 
exists predominantly as the aquated cation. In all cases 
a Single, averaged set of resonances was observed, in- 
dicating exchange between the "free" and complexed forms 
to be fast on the NMR timescale. 

Since the chemical shift of the complexed form 
of methionine could not be eee enaned directly, the for- 
mation constant of the thioether complex was determined 
from the chemical shift data for the methyl protons of 
methionine given in Table V by the iterative procedure 
described previously for the methylmercury valine com- 
plexes, with the exception that all of the "free" methyl- 
mercury at this pH is present as the cation. Using this 
procedure, the formation constant was determined to be 
87 with a standard deviation of 4 and the chemical shift 
of the methyl protons of the thioether-complexed methio- 
nine to be 2.86 ppm. From the chemical shift data for 
the methyl protons of methylmercury and the mercury-proton 
coupling constant data in Table V, the chemical shift of 
the methyl-protons of methylmercury bonded to the thio- 
ether group is calculated to be 1.12 ppm while the mercury- 
proton coupling constant is found to be 223 Hz. 


The formation constant for the amino complex 
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defined by Equations (46) and (47), was determined from 
chemical shift ea ioe the seks proton of methionine 
and the methyl protons of methylmercury for the solution 
containing equimolar amounts of methylmercuric hydroxide 
and methionine. The procedure used was identical to that 
used for valine. Only data at pH values greater than pH 9 
were used to ensure that the only complex formed was the 
amino complex. From the chemical shift data for the methine 
proton, log Kn = 1.0) ¢, 0,02. and the chemical shift of the 
methine proton in the complexed species is 4.21 ppm. The 
chemical shift data for the methyl protons of methylmercury 
yielded a log Kp of 7.40 + 0.08 and a chemical shift of 
0.91 ppm for these protons when methylmercury is bonded 


to the amino group. The individual results are given in 


Tables VI and VII. 


Table VI. Formation Constant Values for Amino Complex 


Calculated from the Shift of the Methylmercury 


Protons 
pH log Kp 
9.00 7.49 
9.50 7.40 
10.00 USS 
10.50 Tosh 
11.00 To43 
7-50 i493 
12.00 Yas: 


3. The Binding of Methylmercury by Cysteine and Penicill- 


amine. 


The binding of methylmercury by the amino acids 


migoion ott 10% 
ab bxonbyd sists it 
sont od La bdieie 
g bie aed? 


eas ape 


ayia io shee : ingle 
oie agg er sotvage 8 | 
yas seins to anos, Says 


ma), 
an B a 


7" to aad octaode. © bie 


‘bebaod, at 


ak wevey. oe ta te 


‘ 
7 aa 
ay 

ty 


zs 


Table VII. 


Formation Constant Values for Amino Complex 


Calculated from the Shift of the Methionine 


Protons. 
pH log Kp 
10.50 Jisoas 45 
10.70 Peo 
11.00 Lipgeks 
fi Ce ons 
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cysteine(III) and penicillamine (IV) was studied to 


HSCH.CHCO., Iit 
+ 
NH, 
HSC (CH,) ,CHCO, LV 
+ 
NH, 


characterize the binding of methylmercury by the sulf- 
hydryl group. The chemical shift of the methyl protons 

of methylmercury in a solution containing 0.150 M methyl- 
mercury and 0.150 M cysteine is shown as a function of 

pH Gin Fig. 13. Also shown for comparison is the shift 

of the methyl protons of methylmercury in a 0.190 M 
methylmercury solution. The mercury-proton spin-spin 
coupling constants for the same solutions are shown in 
Figure 14. The chemical shifts of the resonances for the 
a and £8 carbon atoms of a 0.15 M cysteine solution and 

of a solution containing 0.20 M cysteine and 0.20 M 
methylmercury are shown in igure 15. The displacement 

of the chemical shift and coupling constant titration 
curves of methylmercury when the solution contains cysteine 
and of the carbon-13 chemical shift titration curves of 
cysteine when the solution contains methylmercury indicate 
that methylmercury is complexed by cysteine over the pH 


range 0-14. 
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CHEMICAL SHIFT ,ppom vs TMA 
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Figure 13: pH dependence of the chemical 
shift of the methyl protons of methyl- 

mercury in an aqueous solution containing 
0.190 M methylmercury (solid points) and 
in an aqueous solution containing 0.150 M 
methylmercury and 0.150 M cysteine (open 


poLnts). 
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Figure 14: pH dependence of the mercury-proton 
spin-spin coupling constant of methylmercury in 
an aqueous solution containing 0.190 M methyl- 
mercury (solid points) and in an aqueous solu- 
fiton, containing 0.150 M methylmercury and 0.150 
cysteine (open points). 
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Figure 15: pH dependence of the chemical 
SNL ts Ol tie, =Clos (tA), and «CHN aB)aicarvon 
atoms of cysteine in an aqueous solution 
containing 0.201 M cysteine (solid points) 
and in an aqueous solution containing 
0.200 M cysteine and 0.200 M methylmercury 
(open points). oi 
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To identify which of the three potential binding 
sites is involved in methylmercury complexation, the 
binding of methylmercury by cysteine was studied at pH 
13 since the previous studies with valine apavose ae that 
the amino and carboxyl groups are not methylmercury- 
complexed at this pH. In Figure 16, the chemical shift 
of the methyl protons of methylmercury is given as a 
function of the cysteine to methylmercury mole ratio. 
Also shown for comparison are the results of a similar 
experiment with S-methylcysteine. These results indicate 
that of the 3 potential binding sites of cysteine, the 
sulfhydryl group is the binding site involved in the com- 
plexation of one-coordinate methylmercury at high pH. 
Comparison of the: low pH region of the titration curves 
used in the study of the binding of methylmercury by 
valine (Figures 8 and 9) with the low pH region of Figures 
13 and 14 indicates that the sulfhydryl group is the 
binding site over the pH range 0-13. The mole ratio 
study arSo clearly indicates that only one methylmercury 
is bound by the sulfhydryl group of cysteine at this pH. 

In the “pH ‘region 5-13, the chemical shifts of 
the a and 8 carbon resonances of the 0.150 M cysteine 
solution change as a result of deprotonation of the amino 
and sulfhydryl groups (Figure 15); acid dissociation 
constants for these deprotonation reactions are derived 


from the carbon-13 chemical shift data in Chapter VIII, 
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CHEMICAL SHIFT , ppm vs TMA 
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Figure 16: Chemical shift of the methyl protons 
Of nethylmercury as a “function (ol the rato: of 
cysteine to methylmercury (open points) and the 
ratio of S-methylcysteine to methylmercury (solid 
DOMMES). =<pH = 132.0% 
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Assuming the sulfhydryl group to be deprotonated 
in the cysteine complex of methylmercury, one can at- 
tribute the pH dependence of the chemical shifts of the 
cysteine carbons in the solution containing 0.20 M 
cysteine and 0.20 M methylmercury (Figure 15) to vari- 
ations in the state of protonation of the amino and 
carboxyl groups. Using Equations (41) and (42), the PKay 


for the acid dissociation of the carboxyl group of methyl- 


mercury-complexed cysteine, defined by Equation (48) 


- + 
CHCO,.H + H.O 2 CH.HgSCH.CHCO, + HO (48) 


CH HgSCH 2 2 3 2 2 3 


3 


+ nek 


NH, 3 


was determined to be 1.95 ¥ 0.05 from the chemical shift 
data in the pH region 1-5 and the PK, 5 for the acid dis- 
sociation of the amino group defined by Equation (49) was 
found to be 9.05 t 0.04 from the chemical shift data in 


- - + 
CHCO sae « bm Mead CHHgSCH,CHCO + H30 (49) 


CH 2 2 2 2 


3HgGSCH 


i 
NH, NH, 


the pH range 7-11. The chemical shift of the methyl 
protons of methylmercury in the solution containing 0.20 M 
methylmercury and 0.20 M cysteine (Figure 13) changes with 
pH in those pH ranges where the carboxyl and amino groups 
are being titrated and presumably also reflects the state 


of protonation. The PK > value determined from the data 


80 


———: 
~Kiae Og Tet ome pia i 
ada orikins ort 26. 


st 


bya f shcuiay reo quoty yn oct : is i 


rt am 


Big sot seupe at 9 


’ 


in Figure 13 for the amino group of methylmercury-complexed 


cysteine is 9.0 £.0.1.. The chemical shift data in the 


Al 
of the carboxyl group, possibly because some protonation 


pH range 1-5 did not yield a constant value for the pk 


of the methylmercury-complexed sulfhydryl group occurs 
in this pH range. Protonation of the methylmercury- 
complexed sulfhydryl group is discussed in more detail 
in Section 4 of this chapter. The PK, values for the 
carboxyl and amino deprotonation reactions of methylmercury- 
complexed penicillamine from the chemical shift of the 
methine proton of penicillamine are 2.0 + 0.1 and 9.0 +t 0.1. 
The pH dependence of the chemical shift of the 
methyl protons and of the mercury proton coupling constant 
of the methylmercury in an equimolar solution of methyl- 
mercuric hydroxide and penicillamine is similar to that 
shown in Figures 13 and 14 indicating that complexation 
of the methylmercury to the sulfhydryl group of penicill- 
amine is similar to the complexation with cysteine. These 
two parameters were also measured in basic solutions con- 
taining up to 8 M KOH. Under these extreme conditions, 
the resonance for the methyl protons of methylmercury was 
between that of "free" methylmercury and the methylmercury 
in the penicillamine complex, which may indicate partial 
dissociation of the complex due to competition from 
hydroxide ion for CH3Hg (II) or, possibly «formation of a 


Pate 2- 
complex of the type CHHg (OH) (penicillamine) : 
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The formation constant for the binding of methyl- 
mercury by the sulfhydryl group could not be measured by 


utilizing the competition between H on and CH Hg (IT) for 


3 
ligand or between ligand and OH for CH Hg (IT) Since nowhere 
in the pH region 1-13 does the hydrogen ion compete success- 
fully with methylmercury for the ligand or hydroxide ion 
with the ligand for methylmercury. The determination of 

the formation constant was attempted by using a competing 
ligand which forms a moderately strong complex with methyl- 
mercury. Possibilities suggested by Schwarzenbach and 
Schellenberg's results (2) are thiosulfate ions (log Kp = 


11) and cyanide ions (log K, = 14.1). To pH 6 and 13 


P 
solutions containing equimolar concentrations of methyl- 
mercury and cysteine, thiosulfate was added up to a thio- 
sulfate to methylmercury to cysteine ratio of 4 to 1 tol. 
Even with a large excess of thiosulfate present, no sig- 
nificant amount of the cysteine complex dissociated to 
form the thiosulfate complex, as indicated by the absence 
of a change in the chemical shift of the methyl protons of 
methylmercury. In the cysteine complex, the resonance for 
these protons is at 0.74 ppm whereas in the thiosulfate 
complex it is at 0.97 ppm. 

Similar experiments at a pH of 13 using cyanide 
ion as the competing ligand were somewhat more successful. 


When the solution contained a 20- to 100-fold excess of 


cyanide to methylmercury, the resonance for the methyl 
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protons of methylmercury moved upfield. However, the shift 


of the methyl protons of methylmercury in CH,HgCN differs 


| 
Bye LeSsmecian 20uz, (00S ppm) from the chemical ‘shift in 

the methylmercury-cysteine complex at pH 13. Consequently 
these results could only be used to estimate the logarithm 


of the formation constant for the sulfhydryl complex 


defined by Equations (50) and (51) to be approximately 15. 


- _ + oS - 
SCH,CHCO, + CH,Hg % CH,HgSCH,CHCO., (50) 
NH, NH, 
[CH HoL] 
fae eg ee ae (o1') 
F - + 
[L ] [CHHg ] 


4. The Binding of Methylmercury by Glutathione 

In the previous section it was shown that of the 
potential coordination sites in the amino acids penicill- 
amine and cysteine, the sulfhydryl group is the strongest 
coordination site for methylmercury. Presumably the amino 
and carboxyl groups of sulfhydryl-complexed cysteine and 
penicillamine are still potential coordination sites for 
additional methylmercury. Attempts to characterize the 
binding of methylmercury to these amino and carboxyl groups 
were unsuccessful because the resonances do not indicate 
the state of complexation of a Single functional group. 


To further characterize the binding of methylmercury in 
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such complexes the binding of methylmercury by glutathione- 


(y-L-glutamyl-L-cysteinyl-glycine, V) was studied. 


H O O 
~ O9CCCH.CHCNHCHCNHCH,CO, V 
NH,” CH, 
| 
SH 


The methylene protons of the glycyl residue 


(NHCH,CO. , labelled A in the following discussion), the 


Lie 

methylene protons of the L-cysteinyl residue (CH,-SH, 
labelled B) and the methine proton of the y-L-glutamyl 
residue (0,C (NH,”)CH-, labelled C) reflect the chemical 
state of the carboxylic acid group of the glycyl residue, 
the sulfhydryl group of L-cysteinyl residue, and the 
a-carboxylic acid and a-amino groups of the L-glutamyl 
residue. The chemical shifts of protons A, B, and C are 
shown as a function.of pH for a 0.15 M solution of 
glutathione in Figure 17. The chemical shift of protons A 
changes in the pH range 1-5 one to deprotonation of the 
glycyl carboxyl group, eae of protons B changes in the 

pH range 7-11 due to deprotonation of the sulfhydryl group 
and that of proton C changes in the pH range 1-5 due to 
deprotonation of the a-carboxylic acid group of the y- 


L-glutamyl residue and in the pH range 7-12 due to de- 


protonation of its a-amino group. 
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Figure 17: pH dependence of the chemical 
shattts”’ of “the “(A)“glycyl' methylene protons, 
(B) L-cysteinyl methylene protons, and (C) 


L-glutamyl methine proton of glutathione. 
(0.15 M giutathione, T.= 25°). 
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The chemical shifts of protons A, B and C are 
shown as a function of pH in Figure 18 for a solution 
containing equimolar amounts of glutathione and methyl- 
mercuric hydroxide. A comparison of Figures 17 and 18 
shows that the chemical shifts of protons A and C are 
unaffected by the methylmercury whereas protons B experience 
a marked pH independent downfield shift. This indicates 
that methylmercury is bound by the deprotonated sulfhydryl 
group over the pH range 0-14. | 

The chemical shift of the methyl protons of 
methylmercury and the mercury-proton coupling constant for 
methylmercury in a solution containing equimolar amounts 
of methylmercury and glutathione are shown as a function 
of pH in Figures 19 and 20. Also shown for comparison are 
the chemical shift and mercury-proton coupling constant 
of a 0.190 M methylmercury solution. The data presented 
in Figures 19 and 20 also indicate binding of methylmercury 
by the sulfhydryl group of cutee en ane: 

At pH less than 2, the chemical shift of the 
methyl protons of methylmercury and the mercury proton 
coupling constant shift in the direction of the chemical 
shift and coupling constant of free methylmercury. This 
could be caused by a slight excess of methylmercury or it 
could indicate some dissociation of the complex through 
competition of protons with CH Hg (II) for the sulfhydryl 


group Or some protonation of the methylmercury-complexed 
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Figure 18: pH dependence of the chemical 
shifts of the (A) glycyl methylene protons, 
(B) L-cysteinyl methylene protons, and (C) 
L-glutamyl methine proton of methylmercury- 
complexed glutathione (0.15 M CH3HgOH, 
0.15 M glutathione, T = 25°). 
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Figure 19: pH dependence of the 
chemical shift of the methyl protons 
of methylmercury in an aqueous solu- 
tion containing 0.190 M methylmercury 
(lower curve) and in an aqueous solu- 
tion containing 0.150 M methylmercury 
and 0.150 M glutathione (upper curve). 
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Figure 20: pH dependence of the 
mercury-proton coupling constant of 
methylmercury in an aqueous solution 
containing 0.150 M methylmercury 
(upper curve) and in an aqueous solu- 
tion containing 0.150 M methylmercury 
and 0.150 M glutathione (lower curve). 
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sulfhydryl. The same results were obtained for an ex- 
periment in which the glutathione was present in a 203% 
excess, indicating that the changes in the chemical shift 
of the methyl protons of methylmercury and oe mercury- 
proton coupling constant at pH less than 2 in Figures 19 
and 20 is not due to an excess of methylmercury. To deter- 
mine if these changes are due to dissociation of the 
complex, methionine was added to a pH 0.5 solution con- 
taining 0.17 M methylmercury and 0.17 M glutathione. If 
the methylmercury glutathione complex is partially dis- 
sociated, the methionine would react with the free methyl- 
mercury to form a thioether complex at this pH which would 
be indicated by a downfield change in the exchange-averaged 
chemical shift of the methylmercury protons. No detectable 
change in the methylmercury chemical shift was observed up 
to a methionine concentration of 0.34 M. These results 
suggest that the changes in the chemical shift and coupling 
constant at pH less than 2 are due to some protonation of 
the complexed sulfhydryl group. Such a complex would be 
somewhat similar to the thioether complex formed between 
methylmercury and methionine, in which the mercury-proton 
coupling constant is 223 Hz as compared to a coupling 
constant of 170 Hz in the methylmercury glutathione complex. 
Thus if protonation of the methylmercury complexed sulf- 
hydryl group is occurring at pH less than 2, the exchange- 


averaged mercury-proton coupling constant might be expected 
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to increase as observed. 

Since the S-H stretching vibration in glutathione 
produces a distinctive Raman band in the region of 2580 
cmt, it was thought that evidence for protonation at pH 
less than 2 might be obtained by Raman spectroscopy. Raman 
samples were prepared at pH 0.4 of a solution containing 
0.42 M methylmercury and 0.42 M glutathione. No S-H 
stretching vibration was observed, presumably because the 
concentration of protonated complex is too low to be 
detected by this method. 

To study the binding of additional methylmercury 
by sulfhydryl-complexed glutathione the chemical shifts 
of protons A, B and C were measured as a function of pH 
for a solution containing 0.4 M methylmercury and 0.2 M 
glutathione. Due to overlap of multiplet patterns, it 
was not possible to characterize the binding quantitatively 
with these data. The results suggested, however, that at 
pH less than approximately 7, two methylmercury ions are 
coordinated to the sulfhydryl group. To characterize 
this interaction quantitatively, similar experiments were 
performed by carbon-13 magnetic resonance spectroscopy. 

The carbon-13 magnetic resonance spectrum of a 
0.30 M glutathione solution at pH 7.37 is shown in Figure 
21. Ten well-resolved resonances are observed, one for 
each carbon atom of glutathione. The resonances are iden- 


tified according to the residue and the particular carbon 
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atom within that residue. The chemical shift of each 
carbon atom is sensitive to deprotonation and complexation 
reactions occurring at functional groups in close proximity 
to the particular carbon atom. The chemical shift of the 
individual carbon atoms is shown as a function of pH in 
Figure 22. The resonances for the GLY-C, and GLY-COOH 
carbons ‘shift as the glycyl carboxylic acid group is 
protonated. Deprotonation of the sulfhydryl group is 
reflected in the shifts of the CYS-C,, CYS-C, and CYS-CONH 


carbons. The GLU-C GLU-C., GLU-CONH, and GLU-COOH carbon 


B! 
resonances shift as both the carboxyl and amine groups 
deprotonate. 

To elucidate the nature of higher complexes 
formed from methylmercury-complexed glutathione, the 
chemical shifts of the carbon atoms of glutathione in 
solutions containing methylmercury at methylmercury to 
glutathione ratios of 1:1 and 2:1 were, measured as a 
function of pH. Chemical shift data for the CYS-C, and 


GLU-C, carbon atoms of glutathione are shown in Figure 23% 


B 
Also shown are the shifts of these carbon atoms in a 

solution containing only glutathione. Comparison of the 
chemical shift curve for the CYS-Cp, carbon of free gluta- 
thione (solid points) with that of the CYS-Ce, carbon of a 
1:1 methylmercury glutathione solution (half-open points) 


indicates complexation of the methylmercury by the sulf- 


hydryl group over the pH range studied. The small changes 
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Figure 22: pH dependence of the carbon-13 
chemical shifts of the ten glutathione car- 
bon atoms in an aqueous solution containing 
0.30 M glutathione. 
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in the analogous Gnusc. curves result from changes in the 
acid dissociation constant of the amino group of methyl- 
mercury complexed glutathione when the sulfhydryl group is 
complexed with methylmercury. Comparison of the CYS-C, 
and GLU-C, chemical shift curves for the 2:1 methylmercury 
glutathione solution (open points) with those of the 1:1 
solution indicates that in the pH range 0-4 two methyl- 
mercury cations are bonded to the sulfhydryl group. The 
mercury-proton coupling constant of the methylmercury in 
this species is 212 Hz. As the pH is increased above 4, 
one of the methylmercury cations shifts to the amino group, 
as indicated by the displacement of the GLU-C, 2:1~curve 
from the 1:1 curve. The data in Figure 23 indicate that 
at pH greater than 9, the second methylmercury is partially 
dissociated from the glutathione, presumably to form 
CH3HgOH by analogy with the pH dependence of the binding 
of methylmercury by valine. 

To characterize these complexes further, mole 
ratio studies were performed at pH values 1, 4 and 8. 
The results are presented in Figures 24, 25 and 26. At 
BH] the sbift, of the GLU-C,, carbon remains constant 
(Figure 24) indicating no complexation to the amino group, 
whereas the CYS-C, shift increases up to a mole ratio of 
2 indicating 2 methylmercury cations are coordinated at 
the sulfhydryl site. The break in the CYS-C., curve 


occurs because no averaged resonance could be detected in 
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CHEMICAL SHIFT, ppm vs DIOXANE 


[METHYLMERCURY] / [GLUTATHIONE] 


Figure 24: ..Carbon-13,chemical .shifits..of:the 
GLU-C,, and CYS-C . carbon atoms of glutathione 


as a function of the ratio of methylmercury 
to glutathione. pH = 1.0. 
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CHEMICAL SHIFT, ppm vs DIOXANE 


1 2 
[METHYLMERCURY] / [GLUTATHIONE] 


Figure 25: Carbon-13 chemical shifts of the 


GLU-C., and CYS-C. carbon atoms of glutathione 


as a function of the ratio of methylmercury 
to glutathione. pH = 4.0. 
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CHEMICAL SHIFT, ppm vs DIOXANE 


[METHYL MERCURY] / [GLUTATHIONE] 


Figure 26: Carbon-13 chemical shifts of the 
GLU=C and CYS—Ce carbon atoms of glutathione 


as a function of the ratio of methylmercury 
to glutathione. pH = 8.0. 
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this region due to exchange broadening. The results in 
Figure 25 indicate that also at pH 4 the sulfhydryl group 
is the only coordination site up to a mole ratio of 2. 
However, the GLU-C. chemical shift curve indicates that 
at mole ratios greater than 2 there is some amino co- 
ordination. The results given in Figure 26 indicate that, 
at pH 8, the complexation is somewhat different. Up to 
a mole ratio of 1, there is only sulfhydryl coordination, 
but at mole ratios greater than 1 the second methylmercury 
adds mainly to the amino group as evidenced by the down- 
field shift of the GLU-C. carbon. 
C. Discussion 
1. The Binding of Methylmercury by the Functional Groups 
of Amino Acids and Peptides. 

The NMR results presented in this chapter 
demonstrate that the carboxylic acid, amino, thioether, 
and sulfhydryl groups are all potential binding sites for 
methylmercury with the amount of binding at a particular 
site being strongly pH-dependent. Of these binding sites, 
the sulfhydryl complexes methylmercury most strongly as 
demonstrated by the results in Figure 16 for the mole 
ratio study with cysteine and S-methylcysteine; in the 
methylmercury glutathione complex no dissociation was 
detected over the pH range 0.4 to 13.5. The binding of 


methylmercury by the carboxyl and amino sites is strongly 
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PH dependent due to protonation of the coordination site 

and to formation of the more stable methylmercuric hydroxide. 
Protonation does not compete with complexation of the 
methylmercury at the thioether group and the results in 
Figure 10 indicate that methylmercury binds to the thioether 
group of methionine at low pH where the "free" methyl- 


mercury is present as CH HgOH,”. This indicates that the 


3 
thioether group has a larger affinity for methylmercury 
than does H,0. 

The chemical shift data presented in Figures 23, 
24 and 25 indicate conclusively that the second methyl- 
mercury is complexed by the sulfhydryl group of the 1:1 
methylmercury glutathione complex over the pH range 0.4 to 
6. . The, position of the resonance for the methylene protons 
of the cysteinyl residue of glutathione in a 1:1 methyl- 
mercury guereeni one solution is essentially independent of 
pH, indicating that the methylmercury-complexed sulfhydryl 
group is deprotonated. At pH less than 2, the chemical 
shift of the methyl protons of methylmercury (Figure 19) 
and the mercury-proton coupling constant (Figure 20) 
suggest that there is some protonation of the complex. 

Binding of methylmercury by the sulfhydryl-rich 
protein thionein has been investigated by Chen, Ganther, 
and Hoekstra (26) who concluded on the basis of spectral 


changes, that methylmercury has a lower affinity for 


thionein than does Cd(II). This is opposite to the affinity 
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order obtained from NMR studies of the binding of these 

metal ions by the sulfhydryl group of glutathione (27). 

ame-greater affinity of thionein for Cd(Ii) may result 

from multiple binding of thionein to Cd(II) eenap ane. 

not occur in the Cd(II) complex of glutathione below pH 

6 (27) nor in the methylmercury glutathione complex, 

and presumably not in the methylmercury complex of thionein. 
5 The results in Figure 13 indicate that, upon 

protonation of the amino group of CH3HgSCH.CHNH.,CO, the 

resonance for the methylene carbon which is two bonds 

removed from the site of protonation, shifts by 4.58 ppm 

while that for the methine carbon shifts by only 1.18 ppm 

even though it is one bond closer to the site of protonation. 

Similarly, complexation of the ionized sulfhydryl group 


by CH,Hg(II) causes the methine carbon, which is separated 


3 
from the site of complexation by two bonds to shift by 
2.22 ppm whereas the resonance for the methylene carbon 
shifts by only 0.11 ppm even though it is directly bonded 
to the site of complexation. This is analogous to the 
relative shifts observed in carbon-13 magnetic resonance 
spectra upon protonation of simple amino acids (66) and 


upon complexation of Cd(II) and Zn(II) by the amino group 


of glutathione (27). 


2. The Binding of Methylmercury by Methionine. 


The binding of methylmercury by methionine is 
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of particular interest in view of the previous work on 
the participation of the thioether group in metal com- 
plexes of methionine (63,64,67-69). By a comparison of 
the formation constants of the glycine and methionine 
complexes of Zn(II), Li and Manning (67) showed that 
binding of Zn(II) by methionine is through the amino and 
carboxylate dentates. Ina similar but more extensive 
study, Lenz and Martell (68) concurred with the con- 
clusions of Li and Manning for Zn(II) and proposed that, 
Se sene metal wons) VAg(l):, CalIl), Cd(Z1); CotEL), CutIt), 
Hawi) MG Cll), Mnilt) 7 Ni Lt) ;/Pb(1f1), and Sr tit) only 
Ag(I) binds to the thioether group. McAuliffe, Quagliano, 
and Vallarino (69) showed by infrared spectroscopy that, 
in the solid state at least, binding of several of these 
and other metal ions by anionic methionine is as proposed 
by Lenz and Martell (68). More recently, Natusch and 


Porter (63,64) have demonstrated by proton magnetic 


resonance spectroscopy that a previously undetected complex 


in which Hg(II) is bonded solely to the thioether group, 
forms in acidic solution. Still more recently, Wong, 
Taylor, Chieh and Carty (70) have reported that in the 
complex, D,L-methionatomethylmercury, which was isolated 
from alkaline solution, the methylmercury is bonded to 
the amino group and there is no interaction between the 
methylmercury cation and the thioether group in the solid 


state. 


103 


2 


me. RTO y ayove, at * wed 


Fores) nt hess ond ityioad’ ef 1 aeiaoideon wt d 


oviensine pxom aud te ee fall a Rat 
Meoo: wild aie bexsunage 468) i tte 


spe eth 


pe peerares, bas, (20) boas leramed is 30 


yond iby ats0 ,o% 2 benem sqpom, seveotts 9 os 
tas pieonen fea, texan Am -" ve ie | 


we: he Ok Slat a 
; ip , v 7 : " - Ul 
r Ales | SRA 6 


sepa sedgeodds eae 9) jptetay, ie nto ab ane ro tay mt ‘ 
. é eo Wsaoyer eH08 fina " snodsuloe obbios ab = | ae 
i ads nb aeds besiege: owed (00): we. brie, deh nll a 
povetoat Ba tla rtrosint at : Le 


ena ong : seals | 


Kew (a0 bétosdebaa lene 


=e 


ee oepeos. at _Yausiaemt yea oat: snotsatoa satlexts, mon’ 
oat aioe mo sedihwont: oa a rods, sbemanrne os ts da 


= : 
a a 
ae é is i i s 
‘ 4 t ; ; ' z i Ve 
it } Z fie hw z } ae ps 
: Ll , : SI : 
; a f D Stan its Fj 
v . ae a a fi j eh 
F me { ey 
wil ri # ap | a ble MH! he 
ae : ister : i? iz ' he : a) a) Tr ie. 


The results presented in Figure 10 provide con- 
clusive evidence for methylmercury binding to the thioether 
group.) “The, formation constant for this complex is 
relatively small; thus, it dominates the solution chemistry 
only under VeELV ecidic COondmte1Ons where “the amino 
and carboxylate groups are protonated and are not able to 
compete with the thioether site for the methylmercury 
cations. Because of the high degree of specificity in 
the binding of methylmercury by the thioether group in 
acidic solution, it may be possible to identify methionine 
resonances in the PMR spectra of methionine-containing 
peptides and proteins by observing changes in the spectrum 
as the peptide or protein is titrated with methylmercury 


2VEw (oleh MO par oye 


3. pH Dependence of the Binding of Methylmercury. 

The extent of complexation of methylmercury by 
the molecules discussed in this chapter is strongly pH 
dependent due to protonation of the ligand and reaction of 
methylmercury with hydroxide ion. Thus the conditional 


formation constant, defined by Equations (52) and (53) 


methylmercury, 4, = ligand, 4. 2 Te ee (52) 
Cc [methylmercury oun iex! 
i cr. onan aan ee (53) 


] 


[methylmercury,, J] [ligand,,.. 
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is more informative for a given set of solution conditions 
than formation constants of the type calculated in this 
chapter since the conditional formation constant indicates 
directly the extent of complexation for the particular 

set of solution conditions for which it was derived. The 
conditional formation constant equals aBKoe where a is the 
fraction of "free" ligand in the form which complexes and 
B the fraction of "free" methylmercury as the aquated 
cation for these solution conditions. From the results 
reported in this chapter, the orders of conditional 
stability of methylmercury complexes of sulfhydryl, amine, 
carboxyl, and thioether ligands have been derived as a 
function of pH and are given in Table VIII. HO and OH 
are included to indicate those conditions where the con- 
ditional stability constants are so small that little or 
no complex forms. The information summarized in this 
Table should be useful when considering which type of 
complex will predominate in systems containing multi- 


functional ligands. 


4. Relationship between Formation Constant and Mercury- 


Proton Coupling Constant. 


Ti Chapter’ 1, 1 was noted ‘that, fora series 


of carboxylic acids, the absolute magnitude of the mercury- 


proton coupling constant decreases approximately linearly 


as both the pK, of the carboxylic acid and the logarithm 
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of the formation constant of the methylmercury-carboxylic 
acid complex increase. Evans et. al. (34) and Sheffold 
(71) reported a linear decrease in the magnitude of the 


coupling constant for CH,HgO,CR complexes in deutero- 


2 
chloroform with increase in the aqueous solution PK, of the 
carboxylic acid. In Figure 27 the absolute magnitude of 
the mercury-proton coupling constant is plotted as a 


function of log K,, for thioether, hydroxyl, sulfhydryl 


F 
and amino complexes of methylmercury. The results plotted 
in Figure 27 include those obtained in the present work 
and those obtained in previous studies on these systems. 
An approximately linear relationship is observed. From 
a least squares treatment of the data J = -5.09 log Kp + 


249. This relationship should be useful in determining 


the site of complexation in multifunctional ligands. 
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Figure 27: Plot of mercury-proton spin-spin 


coupling constants versus the logarithm of the 
formation constants of the methylmercury com- 
plex. Values for all ligands studied in this 
thesis as well as values reported in the lit- 
erature are included, ~The’ least squares 
analysis was done using the points indicated 


by arrows. 
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CHAPTER V 


KINETICS OF METHYLMERCURY COMPLEXATION REACTIONS 


Acnakntroduction 

In nearly all of the ligand-methylmercury 
systems presented in Chapter IV, as well as those des- 
cribed in the literature, exchange-averaged spectra were 
observed for both methylmercury and ligand resonances, 
indicating exchange between free and complexed forms to 
be fast on the NMR timescale. This was puzzling, par- 
ticularly in the case of methylmercury complexes of 
sulfhydryl-containing ligands. Since the formation 
constant of a methylmercury complex is equal to the ratio 
of the formation and dissociation rate constants as 


defined by Equations (54) and (55) an upper limit for the 


k 
= = £ . 
Tree CHHg < CH,HgL (54) 
Ka 
k 
c 
K, => (55) 
F ka 


dissociation rate constant can be predicted from the 

formation constant and a diffusion-controlled value for 

the formation rate constant. Such a calculation for a 

methylmercury sulfhydryl complex predicts that Ka will be 
6 -1 


less than approximately 10 sec which corresponds to 


: 6 
a mean lifetime of the complex of approximately 10° sec. 
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The change in the chemical shift of the methylene protons 
of the cysteinyl residue of glutathione upon formation of 
the methylmercury complex is approximately 0.4 ppm so that 
separate resonances for the free and abueiaxea forms of 
glutathione in a solution having a concentration of 
glutathione greater than methylmercury would be expected 

if the lifetime of each form were greater than approximately 
150 sec. Since exchange-averaged spectra are observed, 

the reaction represented by Equation (54) must not be the 
only pathway by which methylmercury and ligand exchange 
between the free and complexed form. In order to elucidate 
the exchange kinetics of methylmercury complexes, the 
exchange kinetics for two systems for which exchange 
broadening was observed were studied. The results are 


presented in this chapter. 


B. Results 
1. The Exchange Kinetics of the Methylmercury Complexes 
of Methylamine. | 

Broadening of the methylmercury resonances has 
been observed in solutions containing amines and some 
simple amino acids over a limited acidic pH range (58). 
In the following section of the thesis, results of a 
study of the exchange kinetics of the methylmercury- 


methylamine complex are presented. 


In a solution containing equimolar concentrations 
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of methylmercury and methylamine above pH 5, exchange- 
averaged spectra are observed for both methylmercury and 
methylamine protons, indicating exchange of both methyl- 
mercury and methylamine between the various free and com- 
plexed forms is rapid on the NMR timescale. Below pH 5, 
the central resonances and the two satellite resonances 
of the methylmercury spectrum are broadened, indicating 
that exchange of methylmercury between the various forms 
occurs at a rate which is measurable on the NMR Riieecale: 
The amount of broadening is pH dependent. In this pH 
region, the resonance for the methyl protons of methyl- 
amine is also broadened; however, it is difficult to 
analyze the broadening in terms of exchange of methylamine 
between free and complexed forms pecauee the methylamine 
proton exchange kinetics are also slow in acidic solution 
and cause broadening of the methylamine resonance. 

The effect of temperature on the NMR spectrum 
of methylmercury is shown in Figure 28. The central 
component of the methyl resonance of methylmercury in a 
solution containing 0.583 M methylamine and 0.146 M 
methylmercury at pH 2.65 is shown in Figure 28. At pH 
2.65, the predominant "free" methylmercury species are 
* and (CHHg) ,0H”. At 1°C, two methylmercury 


CH.HgOH 


is 2 
resonances are observed: the upfield resonance at 0.89 ppm 
versus DSS, arises from the methylamine-complexed methyl- 


mercury, while the downfield resonance at 1.03 ppm is an 
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Figure 28: Central resonance of the methylmercury 
PMR spectrum of a solution containing 0.583 M methyl- 
amine and 0.146 M methylmercury as a function of 


temperature. pH = 2.65. 
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averaged resonance from the "free" methylmercury. Cal- 
culation of the fractional concentrations of methylmercury 
in the free and complexed forms using the reported for- 
mation constant (58) for the methylmercury methylamine 
complex at 25°C indicates that, at this pH, approximately 
one third of the methylmercury is complexed, which is con- 
sistent with the relative intensities of the two resonances 
at 1°C. Thus at 1°C, the exchange of methylmercury between 
free and complexed forms is slow on the NMR timescale, but 
exchange of methylmercury between various "free" forms is 
fast as indicated by the exchange-averaged "free" methyl- 
mercury resonance. As the temperature is increased, the. 
two resonances broaden and coalesce because the rate of 
exchange of methylmercury between free and complexed forms 
increases. At 75°C, a single averaged narrow resonance is 
observed, representing the fast exchange limit. 

In the evaluation of the methylmercury exchange 
kinetics only the exchange broadening of the central methyl- 
mercury resonance was used. Broadening of the satellite 
peaks was also observed but since each satellite constitutes 
only 8% of the total integrated intensity of the methyl- 
mercury resonance, the signal to noise ratio was too low 
to permit the determination of precise lifetimes from the 
widths of these resonances. 

In the methylmercury-methylamine system there 


are several exchange reactions occurring simultaneously. 
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Methylmercury and methyl amine are exchanging between the 
"free" and complexed forms, possibly by the following 


pathways. 


CH Hg! + NH.CH (56) 


+ 
HgNH..CH z 
g cn go5 


3 yeaa) 3 


e a 
CH +ywOH CH,HgOH + NH.CH (57) 


CH.HgNH ve 
aaa 3 yaaa 


3 


+ 


+ 
GHsHGNE CHYs tylb {xt ShCH Hg it. NH.CH (58) 
k 


3 273 3 5 pie: 


+ * * 
HgNH,CH, + NH,CH, = CH,HgNH,CH, + NH,CH (59) 


ied aCHs ee a 7 bape gH, 


3 
In addition, uncomplexed methylmercury is exchanging 

between its various "free" forms and methylamine is ex- 
changing between its protonated and deprotonated forms 


by the reaction 


+ - 
Za 
CH,NH. Fal CH3NH, + H,O (60) 


Of these reactions, only the reactions represented by 
Equations (56) - (58) will cause exchange broadening of 
the methylmercury resonance. Exchange of methylmercury 
between the various "free" forms ‘is fast, as indicated by 


exchange-averaged spectra for methylmercury solutions, 
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while the reaction represented by Equation 59 results 
in exchange of methylmercury between identical environments 
and its rate does not influence the width of the resonance 
line. | 

The lifetime of the complexed species, Tyr was 
calculated from the width of the exchange-averaged central 
resonance for methylmercury over a range of pH values and 
concentrations of methylmercury and methylamine at 25°C. 
The lifetimes were evaluated from the width at half- 
height using Equation (19) presented in Chapter II. To 
evaluate the lifetime from the linewidth for a given 
solution, the populations of the free and complexed forms 
are required, as are the chemical shifts of the central 
resonance for methylmercury in these forms. The chemical 
shift of the complexed form, 2.28 ppm, is known from 
previous work (60). The populations and the chemical shift 
of the free form were calculated as described in Appendix 
I. The results obtained are listed in Table IX. The 
value used for the linewidth of the "free" methylmercury, 
0.38 + 0.04 Hz, was the linewidth of the methyl resonance 
of a solution containing only methylmercury. The value 
used for the linewidth of the complexed species, 0.38 ft 0.04 
Hz was that obtained from a pH 8.5 solution of methyl- 
amine and methylmercury. At this pH, nearly all of the 
methylmercury is in the form of the complex. 


The inverse of the mean lifetime of the complex 
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Table IX. Kinetic Data for Methylmercury-Methylamine 


pH [methyl- [methyl- | OBS ‘i P P 
mercury] amine] m2 oe . f 
M M Hz sec + 

4.51 0.218 O.218 tae 88219 0.282 OC LLS 
4.20 0218 O..218 Les3 84+15 On 2a. OEY ARS) 
3of3 04.218 0.218 £200 Lazo (0 Aa ee 0.849 
3.38 Oe2L8 Oc2b8 1 22 83+16 0.096 0.904 
4.68 Or 218 0.437 L000 80+26 0.543 0.457 
4.23 0.218 0.437 Ly96 60+5 0.404 0.596 
4.01 0.218 0.437 2.02 60+5 05339 0.661 
3.61 0.218 0.437 ZeaU 58+4 0.233 0.767 
B29 ORiaLSo 0.437 1.68 8249 O,163 0.834 
2.80 0.228 0.437 Le bys 116+26 0.085 0.945 
4.10 O75 Q.o2s bes, 69+11 0.456 0.544 
Bide One 0.525 Pe oO A 61+4 0.345 02655 
3.54 O71E75 03525 2.50 60+4 0.275 2.125 
3223" 0.275 0.525 2.06 76+6 0.193 0.807 
2392 OF 175 099525 LAs Lie 1S 02128 O87 2 
4.21 0.146 0.583 122 82+9 0.560 0.440 
390 0.146 OL 58S 1.94 67+6 0.444 05.556 
S768 0.146 02563 2765 5544 OFS 65 0.635 
3.43 0.146 055988 2490 56+4 0.1262 Oe 725 
Deee 0.146 0.583 2.120 1926 O22) VAM STS, 


a Nea 0.146 0.583 1.60 i 225 0.159 0.841 


me eee 


ait .0 cos.0 | esate | Ost at 
efT.0 168.0 pe pote 
epe.0 {2h.0 
$08.0 “880.0 ri 
Ce £2. 9 afede 
ge. ‘BOD rhe 


5 on 
£30.0 — ete.0 inne 208 | ie a _ Feket | 
tae: 9. 665.0 Mae es ae tht 
Meg.0  CaL.g “este” | 8af cones 
218.9 viet aden | ieean ‘aa ‘Petit cn ee 
ee ee ee Pind ie te. e a itr 
829.0. | ane. o- EG oa) 
iia, r ay Sad 
aati. “asta i | aNE0, | 
roe. o eet. me or be Se Ah te ye | 
PSA i ‘wae. reer 7 ff i 

ste.0 98h.0 | eae thet ee | BLO a eS 


ae a é 


Biles 7 
re 


ab. oe.” is ante, om Kale jy be _ fe | 
a22.0 Khe. ee oe Mees £8860 M abo | 2 persia 
ae ages Tin aR Pas, 892.0. ated, ties ” 
eit. SORLm Ae tie “yl $62.0. “ai kes 


'¢ 5 
AED ) ilk » 


erro: “tsa ae aaer les gee Ba. in ‘$8. F pe 
att; 0 eee “ad Boar fd 0 se (fees! 


~—_ 


a8 Aig bo) ae BM - 


zi ree > Ven 
ws oS. pt ¥ ae , it 
¥ i ae ee | ‘s J 


is related to the rate of decrease in the concentration 


of the complex by 


aleltat 
Ae Tel (61) 


where [c] is the concentration of the complex. The rate 
of decrease in the concentration of the complex by the 


reactions which would cause broadening of the methyl- 


mercury resonance [Equations (56) - (58)] 
-d[c) a = a 
rae k, [ce] + k, [OH l[c] + k,[H l[c] (62) 


Division by [c] leads to the following equation which 
relates the experimentally measured mean lifetime to the 


individual rate constants. 


at Y - ; + 
ates k, + k,[0H ] + k,[H ] (63) 


The reaction represented by Equation (58) is 


unlikely to contribute to the exchange since methylmercury- 


complexed amine does not possess a lone pair of electrons 
to promote an electrophilic attack by the proton. The 
reaction represented by Equation (57) can be shown not 
to contribute significantly to the methylmercury exchange 
at the pH values in Table IX by predicting its maximum 
contribution to 1/t, at these pH values. The maximum 
value for ko would be the upper limit for a diffusion- 


10 
controlled bimolecular rate constant at 10 (72) 7 SO 
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that the maximum contribution of the term k, [0H ] at pH 4 
is 1. If the rate constant is less than the diffusion 
controlled value, the contribution would be smaller. The 
magnitude of this term will decrease as vH decreases. 

The uncertainties reported for 1/t, in Table IX were cal- 
culated using only the uncertainty in the observed line- 
widths. Additional sources of error are the values used 
for Pa and Per the fractional populations of free and 
complexed forms but it is difficult to estimate the mag- 
nitude of uthis erxor....This error will lead»to)theelargest 
uncertainty in 1/t, when Po and Pe are very different 
from each other. Since 1/t, does not show, within the 
estimated uncertainties, a systematic increase as the pH 
increases, the reaction represented by Equation (57) must 
not contribute significantly to methylmercury exchange at 
the pH's listed in Table IX. Thus Equation (62) reduces 


to 
4f 


Since the uncertainty in 1/t, due to errors in Po and Pe 
is greatest when Po and Py are very different from each 

other, an average value of 70 2 LS fox ky was calculated 
using only that data where Py is greater than 0.2. From 
the relationship K, = k_y/kyr Ky is calculated to be 

Ee 2 10” bec At pH less than 2, the methylmercury 


resonance is not exchange-broadened since at this pH and 
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lower little or no complexed methylmercury exists. At pH 
greater than 5 the rate of exchange of methylmercury 
becomes fast; presumably because of the contribution of 

a reaction such as that represented by Equation (57) whose 


rate is pH dependent. 


2. The Exchange Kinetics of the Methylmercury Complexes 
of Glutathione. 

With each of the sulfur-containing ligands 
studied, PMR spectra were recorded for the methylmercury 
in solutions having methylmercury to ligand ratios greater 
than one over the pH range 0-13. In every case, a single 
exchange-averaged resonance pattern was observed, in- 
dicating exchange of methylmercury between the various 
methylmercury species to be fast on the NMR timescale. 

However, the rate of exchange of ligand between 
"free" and complexed forms in solutions containing an 
excess of ligand is pH dependent. For the methylmercury- 
penicillamine system, the line width of the exchange- 
averaged resonances Of a solution having a penicillamine 
to methylmercury ratio of 2:1 indicated fast exchange 
except over the pH region 1.7 to 3.5 where a small amount 
of phexdehing of the averaged penicillamine methine proton 
resonance was observed. For the methylmercury-glutathione 
system, the rate of exchange of glutathione between free 


and complexed forms in a solution containing an excess of 
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glutathione is fast at pH greater than 6 as indicated by 
sharp, averaged resonances in both the proton and carbon- 

13 spectra. At pH less than 6, the rate decreases and is 
dependent on pH and glutathione concentration. Due to 

the complexity of the PMR spectrum of glutathione, the 
exchange was studied by carbon-13 NMR. The CYS-C region 
Pfeene Carbon-13 spectra, for a solution containing 0.33 M 
Siutachione jand 0.165 Mimethylmercury is shown as a function 
SrepHy in, Pig.7295) the eo. resonance is shifted more 

upon complex formation and thus is more sensitive to kinetic 
effects, however overlap of the exchange-averaged CYS-C 


B 


resonance with the GLU-C, resonance precludes analysis of 
its lineshape. At these concentrations, the exchange rate 
Demsnrtrcrentiy slow atypH 2.12 that, separate exchange-— 
broadened resonances are observed. 

The inverse of the mean lifetime of the methyl- 
mercury-glutathione complex is given in Table X for a range 
of experimental conditions. The mean lifetimes were 
obtained by matching experimental and computer-simulated 
spectra. Spectra were simulated as a function of the 
lifetime of glutathione in each of the environments using 
the Bloch phenomenological equations as discussed in 
Chapter ible milhe wuncentainty sofi the (lifetimes; ;duewmainly 
to the low signal-to-noise ratio for the exchange-broadened 


Spectra, of which thoselan Figure 29 are typical, is 


estimated to be ~t 153%. 
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Figure 29: CYS-C, region of the carbon-13 magnetic 


resonance spectra of a solution containing 0.165 M 
methylmercury-glutathione and 0.165 M glutathione 
as a fuMmeiionsor pH. -25-C. 


oa 


) Bese ti-nediss edd to “hobpex ‘anes J-aND 48 pauls ne | 
ae He ‘paiatszxop: nod vies 6 to ei siekaoee oe 

ir | Pi Lieeegcogn a 

| i 


122 


Table X. NMR Kinetic Data for the Methylmercury-glutathione 


System = 
pH [complexed [free Be = 
glutathione] glutathione] Cs 
M M sec | 
2493 Vero 0.165 100 
Be52Z OsLG5 0.165 33 
22352 0.165 0.165 25 
ealZ O65 05165 17 
reoo 02165 0.165 20 
D5 0.165 0.165 a3 
16 0.165 O-165 50 
ae0S On165 0.165 40 
0590 On 165 0.165 50 
OTS OeLe5 Oi L65. = 200 
0265 ») Oal65 0.165 200 
0.48 0.165 O7.1655 240 
ZLo8 0.28 0.28 67 
2.54 0.28 | 0.28 50 
Case 0.28 0.28 50 
2.14 0.28 O28 40 
Maes: 0.28 0.28 33 
1.45 0.28 O.28 40 
2.06 0.28 0.28 50 
095 0.28 | 0.28 50 


(Table continued) 
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Table X continued 


pH [complexed 
glutathione] 

M 

0.80 0.28 

0.68 0.28 

0.47 0228 

Gr2Ssc 

b 


simulated spectra. Uncertainty in 1/t, values is estimated 


to. be + 153. 


[free 
glutathione] 


200 


Obtained by comparison of experimental and computer 
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Possible reactions by which glutathione might 
exchange between the free and complexed forms for the 


conditions of Table X are: 


1 
CH,HgSG 2 CHjHg’ + GS” (65) 
-1 
k 
+ ot + 
CH,HgSG + HO 2 CH,Hg + GSH (66) 
k 
-2 
CH,HgSG + OH < CH,HgOH + GS (67) 
k 
-3 
k 
ae a * = 
CH,HgSG + GS re CHHgSG + GS (68) 
k 
* 3 * 
CH,HgSG + GSH < CH,HgSG + GSH (69) 


where GSH, GS, and CH,HgSG represent sulfhydryl-protonated, 
sulhydryl-deprotonated and sulfhydryl-complexed glutathione. 
The rate of decrease in the concentration of the complex 

is given by Equation (70) 


A[C} 
dt 


+ - 
= k, [C] + k, [H }{c] + k, [OH LEC} + k,¢o IF] [Cc] + 


kpog,,[F) (C] (70) 


where [F] and [C] are the concentrations of free and 


complexed glutathione and ao = [GS ]/{(F] ana Aon = [GSH] /[F]. 
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Division by [C] leads to the following equation for the 


inverse of the mean lifetime of the complex. 


ies a = 
ap =k, + kj [H ] + k,[0H ] + k,o,[F] + k Eee ol) 


5(sH 


The relative importance of exchange by the reaction 
represented by Equation (65) can be predicted. By 


assuming that the reaction of CH Hg* and GS in Equation 


3 
(65) is diffusion-controlled and therefore the rate con- 


1 has a value of toe° Me cece we predict that 
ky = k_jK, = nies Salo: using the formation constant 


reported by Simpson (20) for the methylmercury-glutathione 


stant k_ 


complex. If the rate constant Ky is less than diffusion- 


controlled, k, will be even less, and thus the reaction 


1 
represented by Equation (65) is predicted to contribute 
a negligible amount to 1/t.- The reaction represented 
by Equation (67) can also be shown not to contribute 
significantly to the observed exchange, at least up to a 
pH of approximately 4. If the reaction of OH with 
HgSG is diffusion-controlled and has a rate eonstant 


3 
LOT = i ot ee : : 
of 10 Ma sec, Les, contribution to 1/t, would 


CH 


only be 0.01 sec + at pH 2 and l sec} at pH 4. It is 


likely, however, that k, is somewhat less than diffusion- 


4 
controlled since Simpson reported a value of 1.6 x 10 


mM tsec! for the rate constant for the displacement of CN” 


from CHHgCH by OH (32) in which case the contribution 
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to V/t, from the reaction represented by Equation (67) 
will be even less. 

The relative importance of the reactions rep- 
resented by Equations (66), (68) and (69) can be predicted 
from their dependence on solution conditions. The spectra 
shown in Figure 29 and the more extensive data in Table X 
for the solution containing 0.165 M methylmercury- 
gluthathione and 0.165 M free glutathione indicate that 
the rate of exchange is pH dependent; the minimum rate 
being at pH 2.1. The larger exchange rate at pH values 
less than 2.1 indicates that exchange occurs by proton- 
assisted dissociation of the complex at these pH values, 
while the increase in exchange rate at pH greater than 2.1 
is due to displacement of complexed glutathione by Gsisr 
whose concentration increases as the pH increases. If the 
rate of exchange by the pH dependent reactions represented 
by Equations (66) and (68) is negligible at pH 2.12, the 
value of 1/t, at this pH is equal to keto fF]. The rate 
of exchange by the reaction represented by Equation (69) 
will be pH independent for che pH range of the data in 


Table X since a the fraction of glutathione present as 


SH’ 
GSH, is greater than 0.999 over the pH range 1-6. There- 
fore the contribution of the reaction represented by 
Equation (69) will be the same as that at pH 2.12. If, 


however, some exchange is also occurring by the reactions 


represented by Equations (66) and (68) at pH 2.12, the 
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contribution of the reaction represented by Equation (69) 
to the observed 1/t will be even less, indicating that 
the amount of exchange by the reaction represented by 
Equation (69) for the conditions given in Table X is 
likely to be small relative to that by Equations (66) 

and (68). Since this value of 17 sééek for V/t, at pH 


2.12 represents the upper limit for k spl FI and is small 


5 

compared to the variation observed in V/t or we shall 

neglect the contribution of this reaction to 1/t, ina 

preliminary analysis of the data. Equation (71) then 

reduces to 

re i |p kde LE] (72) 
2 47S 


T 
Cc 


The exchange-broadened resonances of a solution containing 
0.56 M glutathione and 0.28 M methylmercury never separate 
into the resonances of the free and complexed forms, in- 
dicating that the exchange rate is concentration dependent. 
The concentration dependence is also indicated by line 
shape changes Meneie ratio experiments. 


4 


data in Table X by a method of successive approximations. 


Rate constants ko and k, were estimated from the 


First, a value was calculated for k, from each of the 
lifetimes at pH less than 2 by using Equation (72) and 
assuming k5[H"] >> k,a,(F]. This procedure yielded an 
average value of 1.27 x 102 fOr k,- Then a value for ky 


was calculated from each of the lifetimes for pH greater 
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than 2.1 using Equation (72) and this average value for 
k,- Table XI lists the results of calculations of a new 


value for ky using the mean value for k The mean value 


Ae 
of k, is found to be 600 Meesecu- with a standard 


deviation of 200; this value was then used to recalculate 


k,- The mean value of the results listed in Table XII 


8 a ib 


for k, is 5.8 x 10 M ~sec with a standard deviation of 


4 


1? Seg 10°. From the value of ko, k_» is calculated as 
9 


De. x LO m ‘sec + using the relationship 


PS rd ws (73) 
where Koy is the microscopic acid dissociation constant 
(fo). Lor the sulthydry) group and Kp is the formation 


constant (20). 

In the calculation of the rate constants it was 
assumed that exchange Sse the reaction represented by 
Equation (69) was negligible. If however, the exchange 
at pH 2.12 is due solely to the reaction represented by 


Equation (69), the rate constants obtained by the above 


procedure would be k, = 400 + 250 Ma secune ee 
ere et0  secice ee liane iG Ghat sacs ana 
k. “=. 110 Mm sectt. Since the actual contribution of the 


5 
reaction represented by Equation (69) probably lies some~ 


where between the two extremes used in the above cal- 
culations, the two sets of rate constants can be considered 


as upper and lower limits on the actual rate constants. 
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Table XI. 


Calculated Values fork... 


Expt. 


mole ratio 


0..3/0.15 
0.3/0 %h5 
O33 / Ole Lo 
02370..55 
0238/0 eo 
053/0.15 
025/0.20 
W/O 20 
625/0...25 
0570225 


O25/70..20 


pH 


2 
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Table XII .<«. Calculated. Values. for. k,. 


4 
ely k 
Experiment pH 7 [F] c 4 
123 é A Pe 
sec7l M sec 
0.3/0.15 3.42 roar’ Pea 100 rae rele 
i370. 15 2.93 es stOT 100 65100 
0,370.15 2.52 5.88x10°° thee 5.4x10° 
0.3/0.15 2.32 cia eo me 25 5.9x10° 
0370.15 2.12 B Oexl07” 16.7 Silo 
0.5/0.25 2.93 eens 66.7 2.6x10° 
0.5/0.25 2.54 103x108 50 4.7x10° 
0.5/0.25 2.32 6 2l0ne 50 Teco 
0.5/0.25 2.14 be Seawap 40 8: 8x10° 
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C. Discussion 

For most methylmercury systems, the observed 
NMR spectra indicate that ligand and methylmercury ex- 
change is fast on the NMR timescale. As mentioned in 
the introduction to this chapter, this was somewhat sur- 
prising considering the stability of some of the complexes 
formed. These observations however can be explained in 
terms of two aspects of the complexation chemistry of 
methylmercury as demonstrated by the kinetic results for 
the two systems studied in this chapter. 

Firstly, methylmercury has a tendency to form 
polynuclear complexes. Reactions of the following type 


occur 


CH,HgX + CH Hg" (CHjHg) 9X" (74) 


3 f 


which have the effect of labilizing the CH,Hg-X bond. 

In the methylmercury-methylamine system, 
exchange was slow enough to be observable on the NMR 
timescale because complexes of the type (CH,Hg) 5X cannot 
be formed. However, as shown in the previous chapter, a 
second methylmercury can bind to a methylmercury-complexed 
sulfhydryl group which apparently labilizes the sulf- 
hydryl-complexed methylmercury and results in fast, 
exchange-averaged spectra for Soiibedionee having a methyl- 
mercury to glutathione ratio greater than l. 


Secondly, methylmercury complexes can react 
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with excess ligand according to the following reaction 


CH ,HgX ee fag CH ,HgX* + X (75) 


These two types of reactions have the effect of causing 
exchange to be very rapid whenever there is even a slight 
excess of methylmercury or of ligand. 

In the case of methylmercury-glutathione, fast 
exchange is observed for both proton and carbon-13 
resonances except when excess glutathione is present. 
‘Results presented in this chapter show that reactions of 
this ligand displacement type are the predominant pathway 
for exchange of glutathione between free and complexed 
forms. Above pH 6, exchange is rapid. Below this pH, 
the concentration of sulfhydryl deprotonated ligand is 
small and the exchange by this pathway is less. At pH 
less than 2, exchange is rapid again due to a proton- 
assisted dissociation of the complex. In Chapter IV, 
results were presented which suggest the existence of a 
protonated-methylmercury-complexed sulfhydryl which 
possibly is an intermediate in the proton-assisted dis- 
sociation of the complex. 

These conclusions are significantly different 
from those of Simpson, Hopkins and Hague (28) who studied 
the related ligand N-acetyl-L-cysteine. They concluded 
from PMR measurements of methylmercury-N-acetyl-L-cysteine 
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forms is by a pathway involving first order dissociation 

of the methylmercury complex. Their conclusions are 

questionable however, since the rate constant for the 

formation of the complex is predicted from their dis-— 

sociation rate data to be six to seven orders of magnitude 

larger than diffusion—controlled bimolecular rate constants. 
Bigen, Geier and Kruse (4) have studied ligand 


displacement reactions of the type 


k 
CH,HgOH + ae 2 CH Hgx + OH (76) 
k 


where X is a halide ion. They have shown that E55 is 


strongly dependent on the identity of ZX, while Koa is 
reasonably independent of ZX and they have proposed that 
the reaction proceeds by an associative mechanism in- 


volving an intermediate of the type 


The formation constant of CHLEGO= is larger than those 
of the halide complexes, suggesting that, once such an 
intermediate forms, the tendency for X to dissociate 

from the intermediate will be larger than for OH . In 
support of this (2) E5> < Ea) for the halides (b) both 


K, and k,, increase in the order Cl < Br <I, and 
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(c) when X is CN. which forms a methylmercury complex 


having log Ka = 14.1, compared to log Kp = Ons Thor 


3HgOH, Lal is less than Kio and is approximately 5 


orders of magnitude less than when X is a halide ion. 


CH 


The mechanism of reactions of the type rep- 
resented by Equation (75) presumably proceeds through an 
intermediate analogous to that proposed by Eigen, Geier 
and Kruse for methylmercuric hydroxide reactions. If 
this is the case, there is an equal probability that once 
the intermediate forms, it will dissociate to products or 
reactants. Thus the rate of ligand exchange is predicted 
to be one half the rate of formation of the intermediate, 
which is expected to be determined by the rate of dif- 
fusional encounter or to be slightly less than diffusion- 
controlled rate constants due to steric effects. 

It is not known if complex formation reactions 
involving CHjHgOH,” proceed by an associative or dis- 
sociative mechanism. In the formation of the methylmercury- 


methylamine complex from CH ,HgOH,” and CH.,NH which has 


Sha epee! 
a rate constant of 3 x is Bear the possible mechanisms 
are: for an associative mechanism, the reaction would 


proceed through an intermediate such as 
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While for a dissociative mechanism the reaction would 


proceed through the following two steps 


CH HgOH,” 2 cH 


+ 
H ep C7t) 
3 g H,O 


3 


x hey 7 78 
CHjHg + CH,NH, 2 CH,HgNH,CH, (78) 


with the first being the rate-determining step. The 
results presented in this thesis for the reaction of 
methylamine with CHHgOH,” and other results for CH,HgOH,” 
reported in the literature do not allow one to distinguish 
between these 2 possible mechanisms. With other metal 
ions, the mechanism has been identified by studying the 
dependence of the rate of exchange on the identity of 

the ligand. For example, if the rate of a complexation 
reaction for an aquated metal ion M(H,0) does not 
depend on X, the mechanism is assumed to be dissociative 
with the rate of water loss being the rate-determining 
step. On the other hand, if the rate depends on X, an 
associative mechanism is assumed. The data presently 
available for complexation reactions of CHHgOH.,” do not 

uy reactions because, 


allow such an analysis for the CH,HgOH 


3 2 
for the systems studied to date, the formation constants of 
the CH, 


are associative and proceed through intermediates such as 


CH Hg? ~ , the rate of the reaction will not depend 


HgX complexes are so large that even if the reactions 
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on the nature of X since dissociation of water from the 
intermediate will be faster than dissociation of X. One 
way to resolve this problem might be to study ligands 
whose affinity for the methylmercury cation is similar 
to that of water. Ligands that might be suitable are 
those with the thioether functional group. 

The results reported in this chapter for the 
ligand exchange kinetics in the methylmercury-glutathione 
system have implications as to the lability of methyl- 
mercury in biological systems. In biological systems, 
methylmercury is bound by sulfhydryl-containing proteins. 
In the sulfhydryl-containing ligands studied in this 
thesis the methylmercury binding was very labile. This 
was shown to be due to the presence of excess methyl- 
mercury or excess sulfhydryl ligand. When methylmercury 
is bound by a protein it is unlikely that excess methyl- 
mercury is present or that there would be another nearby 
sulfhydryl group in a position suitable for binding. 

Thus these results suggest that the methylmercury bound 


by a protein is probably quite inert kinetically. 
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APPENDIX I 
Calculation of Concentrations of Species Present 


in Methylmercury-Methylamine System 


The following derivation is for the calculation 
of the concentrations of the species present in solutions 
having equal total concentrations of methylmercury and 


methylamine, i.e., 


(CHHG) moran, = (MeAm) mopar arena a 
The formation constant, Kar of the methylamine complex is 
known from previous work (58): log Kp at le 

K be [complex] (2) 

F [CH,NH,] [CH Hg (1I)] 


Letting gs equal the fraction of methylmercury which is 
complexed, Equation (3) follows from Equations (1) and 


(2). 


(PL) 
Ke! = Qagueegis ears: (3) 
B ew? On eae 
[CH.NH. ] 
where a, = eee (4) 
[CH,NH, ] + [CH,NH, J 


Expressing the concentrations in terms of the acid dis- 


sociation constant leads to 
a = — (5) 
[H ] + Ky 


In Equation (3), 8 represents the fraction of methyl- 
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mercury in the cation form. 

The free methylmercury is distributed between 
three forms: CH,Hg (II), CH, 
trimer, (CH3Hg) ,0", is ignored since calculations show 


HgOH, and (CH3Hg) OH”. The 


that only a very small amount is present at these concen- 
trations of methylmercury (see Chapter III). Because of 
the unsymmetrical nature of the equilibria relating the 
concentrations of the three forms of methylmercury, the 
concentration of methylmercury in the cation form depends 
on the concentration of free methylmercury. Thus Fe and 
B must be determined simultaneously since the value for 
one depends on the value obtained for the other. 


Equation (3) may be rewritten as 


F oUp 
pe sore ee Tt 7B c rr “Bp 


which rearranges to 


2 2 o 
= Op 8 + eg + KeLp a, 8 =O) (7) 


Z 


2 
(Koop a,8) P - (2K,L 


Equation (7) is a quadratic of. the form ax? + bx +c = 0 


2 
where a = K LL 1,8 
Ba oe (2 aT ae B. +. L-) 
Fr DPasB a 
2 
c = KoLp 0.48 
Thus 2 
-b 2 b i 4ac (8) 
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Only the negative root of Equation (7) is physically 
meaningful. The method of solving for P, was to make an 
initial guess at 8, the fraction of methylmercury in the 
cation form. A suitable initial guess for ® is 8 = 1 which 
is equivalent to assuming that all the methylmercury is 

an, the cation, form. oe is then calculated using Equation 
(8). Then using a concentration of free methylmercury 


equal to (1 - PAL a new value of 8 was obtained using 


-! 
Equation (38) Chapter IV. This new value of 8 was then 
used to calculate a new value for Por The iterative pro- 
cedure was repeated until Pa and 8 both converged. 

Once P. and 8 are known the concentrations of 


the various forms of methylmercury can be calculated 


using the following equations: 


a ae = 8 (1-P Lp (9) 
~4.70 - 
[CHHgOH] = 10 ({[CH,Hg (II) ]) ([OH }]) (10) 
21 (CHyHg) OH] = L» - [complex] - [CH,Hg(II)] - [CH,HgOH] 
Cha) 


From the concentrations of the three forms of free methyl- 
mercury and their chemical shifts, the exchange averaged 
shift of the free methylmercury was calculated using 


Equation (12) 


= + 6 
5 aa PoATION CATION a P DIMER SDIMER PBASE BASE 
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where 


[CHHg (II) ] 


PCATION = (13) 


[CH,Hg(II)] + [CH,HgOH] + 21 (CHjHg) OH") 


21 (CH3Hg) OH] 


[CH,Hg(II)] + [(CH,HgOH] + 2[ (CH3Hg) , OH ] 


[CHHgOH] 


P 


BASE (15) 


[CH,Hg(II)] + [CH,HgOH] + 2[ (CH Hg) OH") 


Using y and Sar the lifetime of the complexed species 


was calculated as described in Chapter V. 
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PART 2 


THE ACID-BASE CHEMISTRY OF AMINO ACIDS 
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CHAPTER VI 


INTRODUCTION 


The acid-base chemistry of amino acids has been 
a subject of interest for some time; acid dissociation 
constants have been used in the identification of amino 
acids and in characterizing their solution chemistry. 

The acid-base chemistry of amino acids is 
quantitatively described by stepwise acid dissociation 
constants; for example, for a triprotic amino acid such 


as fully-protonated cysteine, 


HSCH., CHCOOH 


+ 
NH, 


HAL’, acid dissociation occurs in the following stepwise 


fashion and is described by acid dissociation constants 


Ky, K, and K.- 
+ + 
H5L z H + H5L (1) 
[H"] [HD] 
Ky = PRIS (2) 
[HL ] 
> Ay a 
HL a +O (3) 
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[H*) [HL] 
K, = ——_—_ (4) 
[H,0] 
BECH 2H i" (5) 
fy) (a-") 
Ky = (6) 


[HL ] 


The most common method for the evaluation of 
acid dissociation constants is based on pH titrations. 


By comparison of the K, values of cysteine with those of 


A 
carboxylic acids, amines and sulfhydryl compounds it has 


been established that K, describes the deprotonation of 


1 
the carboxylic acid group. “The assignment of K, and K3 
to the remaining acidic groups, however, is not so 
straightforward. The PK, for the sulfhydryl group of 
ethylmercaptan is reported to be 10.50 (1) while that for 
the amino group of methylamine is 10.72 (2). Thus the 
acid strength of the ammonium group is of the same order 
of magnitude as that of the sulfhydryl GLOUD. Con 
sequently, at the molecular level, acid dissociation of 
deprotonated cysteine can occur by two different path- 


ways, which are shown below in the microscopic acid 


dissociation scheme. 
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i. NH,” K193 
HS tH: e if x 
“ee — ee 2 Hl ae 
+ + 
NH, 


NH, ye NH, 
(1) as Eis anuiay) 
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NH, 


(ELT) 


The microscopic acid dissociation constants, defined by 
Equations (7) - (10), are generally represented by lower 
case k's (3). The acid dissociation step which is des- 
cribed by a given microscopic constant is specified first 
by the last number in the subscript, which denotes the 
position of the proton involved in the reaction under 
consideration, and second by the whole set of preceding 
numbers in the subscript, which, regardless of order, 
denote the other positions from which protons have already 
been removed. In cysteine, the carboxyl, sulfhydryl and 


ammonium groups are denoted by numbers 1, 2 and 3 res- 


pectively. 
+ 
=. Pai 
Kp Giwey ir] i 
+ 
et ber] (8) 


ky3 = [T] 


5 : + 
“a = a an i 
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Niee a) TrrtT (10) 


It is not possible to determine the microscopic 
acid dissociation constants of cysteine from its pH tit- 
ration curve alone. By pH titration the concentration 
of the monoprotonated form, which is the sum of the 
concentrations of the two individual monoprotonated forms, 
can be determined but not the concentrations of the two 
monoprotonated isomers. 

The macroscopic acid dissociation constants 


can be shown to be composites of the microscopic constants, 


el (HOepET) eecraTEY 


ae (I) Oe 
Bore io) wena 2) 
*3 ARG OR es 


3 123 132 


Several pathways are available for the de- 
protonation of several of the common amino acids, A few 


examples are lysine, aspartic acid, glutamic acid, 
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penicillamine, tyrosine, and homocysteine. Many peptides 
also contain several functional groups of similar acid 
strength and thus can be characterized in this way. For 
example, any peptide containing a lysine residue has at 
least two amino groups, the N-terminal amino group and 
the lysine-6-amino group. Peptides containing glutamic 
or aspartic acid residues contain a terminal carboxylic 
acid plus the carboxylic acid group in the glutamyl or 
aspartyl side chain. Glutathione is an example of such 

a peptide. A knowledge of the microscopic dissociation 
constants for such peptides may be useful when studying 
the chemistry of these molecules. For example, in con- 
Sidering the coordination chemistry of a molecule such as 
cysteine, if one of the two monoprotonated forms spec- 
ifically forms a metal complex, the microscopic constants 
must be known in order to determine the formation constant 
of the complex. 

In general, most compounds have not been 
characterized in terms of their microscopic dissociation 
constants because the required information is not avail- 
able from pH titrations. To determine the microscopic 
dissociation constants a method of monitoring dissociation 
at specific sites is needed. The various methods which 
have been used to determine microscopic ionization con- 
stants will now be described using cysteine as an example 


since it has been studied the most extensively. 
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In oie oes system, a popular approach has 
been the use of model compounds. Ryklan and Schmidt (4) 
assumed that the effect of an Sealkyl group on the de- 
protonation of a neighboring group would be approximately 
the same as that of an un-ionized SH group, They deter- 
mined the PK, for the ammonium group in S-ethylcysteine 
to be 8.60 at 25°, which they considered to be equal to 


the microscopic constant pk for cysteine itself. They 


13 
then used this value in conjunction with the macroscopic 
constants to determine values for the remaining micro- 
Scopic constants. . [Equations @b2)siandm(d4)i)ia Their 
results are listed in Table XIII. 

Grafius and Neilands (5) determined the micro- 
scopic dissociation constants in a similar way using the 
PK, determined for the sulfhydryl group of cysteine 
betaine, which contains a -N(CH3) .” group instead of the 


-NH i group of cysteine. These workers assumed the 


3 
trimethyl ammonium group resembles the NH,” group in ues 
effect on the deprotonation of the sulfhydryl group and 


of cysteine betaine equals pk of 


thus, thatthe pk 12 


A2 
cysteine. 

A more direct approach to the problem was made 
by Benesch and Benesch (6), who used the ultraviolet 
absorption spectrum of the ionized -SH group. From the 


ultraviolet data, the fraction of the cysteine sulf- 


hydryl groups in the deprotonated form could be determined 
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Table XIII. 
Pki5  PKi3 
S266 8.60 
B65 += 8.75 
eySayl 8.86 
o950 11 -6..85 
Bis Sear 
Geb 4)» 8.86 


Cysteine. 


PK 93 


10.45 


LO: Oi 


PKi 39 


10.51 


method 


pH titration using 
S-ethylcysteine as 
model compound 


pH titration using 
S-methylcysteine and 
cysteine betaine as 
model compounds 


ultraviolet absorp- 
tion spectroscopy 


Raman spectroscopy 
calorimetric methods 


ultraviolet spect- 
roscopy 


Microscopic Acid Dissociation Constants of 


ref. 


(4) 


(5) 


(6) 


(7) 
(8) 
(9) 
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from the ratio of the absorption at the given pH value 

to the absorption when the -SH group is completely dis- 
sociated. Using this fraction along with the macroscopic 
dissociation constants, the microscopic constants were 
g@ecermined. "Their results are listed in Table XIIT, 

The acid-base chemistry of cysteine has been studied by 

a number of other workers using the ultraviolet absorption 
method developed by Benesch and Benesch. Some of their 
results are also listed in Table XIII. 

Elson and Edsall (7) determined the microscopic 
constants by monitoring the deprotonation of the sulf- 
hydryl group by Raman spectroscopy. The results they 
obtained at relatively high ionic strengths are close to 
those obtained by Benesch and Benesch who worked at an 
ionic strength mear *O%.2% 

Wrathall, Izatt and Christensen (8) described 
a calorimetric method for determining microscopic acid 
dissociation constants. The method is based on the fact 
that the AH values for proton dissociation from the NH,” 
group of eae isicysten we and the -SH group of mercapto- 
acetic acid differ by approximately 4 kcal/mole. The 
values they reported varied greatly with ionic strength; 
their results at an ionic strength of 1 are reported in 
Table XIII. This method has been criticized by Wilson 
and Martin (10), who point out that the calorimetric 


method requires that heats for both sulfhydryl de- 
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protonations and, separately, both ammonium deprotonations 
be equal and hence independent of the ionic state of the 
remainder of the molecule, 

| Most of the experiments described ‘above for 
the determination of the microscopic acid dissociation 
constants of cysteine were possible because of the 
special properties of the sulfhydryl group. Thus these 
methods are not applicable to many amino acids and pep- 
tides, for example amino acids or peptides containing 
two amino groups or two carboxyl groups, There is no 
direct method available for systems of these types, 

Since the chemical shift, 64, of a magnetically- 
active nucleus is sensitive to the state of protonation 
of nearby acidic groups, nuclear magnetic resonance 
spectroscopy might be a general method for determining 
microscopic acid aiesectavion constants. Grunwald, 
Loewenstein and Meiboom (11) found that, in solutions of 
methylamine and methylamine hydrochloride, the chemical 
Shift of the methyl protons 1s"a’ linear Tunction’ of’ the 
concentration ratio of acid to base plus acid and thus, 
that it is possible to Seen acid dissociation constants 
by using the chemical shift to determine the fractions 
of protonated and deprotonated forms. 

Proton magnetic resonance spectroscopy has been 
used by several groups of workers to measure microscopic 


dissociation constants. Lowenstein and Roberts (12) 
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studied the ionization of citric acid using the shift of 
the methylene protons as a measure of the state of 
protonation. To determine the effect of the central 
carboxyl deprotonation on the methylene shift they used 
the citric acid symmetrical monomethyl ester. Thus they 
assumed the methylated carboxyl group behaves with regard 
to the chemical shift as if it were a non-ionized carboxyl 
group. Their results have been criticized by Martin (13). 

Rigler et. al. (14) investigated the microscopic 
acid-base equilibria of Boar belies using proton magnetic 
resonance techniques. They used a quaternary methyl- 
ammonium derivative to assist in separating the effects 
of deprotonation of various grups on the chemical shift. 

More recently, Walters and Leyden (15) have 
applied proton Taameeae resonance spectroscopy to the 
problem of microscopic constant determination of cysteine. 
Since the shift of the protons is affected by both amino 
and sulfhydryl deprotonations they found it necessary to 
use a model compound, S-methylcysteine, to assign shifts 
to the intermediate monoprotonated species. The validity 
of such a procedure is questionable and is discussed 
furtherjin Chapter VIELE) otheir aesults careim close 
agreement with the calorimetric results (8) but disagree 
with the results obtained by other methods (6,7). 

In all examples of the application of NMR to 


microscopic constant determination discussed thus far, 
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the chemical shift was affected by more than one de- 
protonation reaction occurring simultaneously. Model 
compounds were used in each case to separate the effects 
of these deprotonation reactions so that chemical shift 
values could be obtained for the intermediate species. 
The situation is simplified somewhat when the chemical 
shift is affected only by deprotonation at a single site. 
There are a few examples of such studies in the literature. 
The acid-base chemistry of fully protonated glutathione 
and methylmercury-complexed glutathione has been studied 
by PMR (16). ° In glutathione, the simultaneously- 
deprotonating groups are well removed from each other, 
thus the chemical shift of a given proton reflects the 
fractional deprotonation of only the nearby acidic group. 
The microscopic constants were determined from the 
fractional deprotonation of each of the acidic groups 
asa function, of: pH. 

Creyf,-.van Poucke , andsBeckhaut: (17) used, a 
similar method to obtain the microscopic acidity con- 
stants of deprotonated asymmetric N-methylsubstituted 
ethylenediamines and N-methylpiperazine. The shift of 
the methyl group is influenced only by changes in the 
protonation of the nitrogen atom adjacent to the methyl 
group. 

These workers all used proton magnetic res- 


onance techniques. Carbon-13 magnetic resonance is a 
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potentially more useful tool, In the proton-decoupled 
carbon-13 spectrum single lines are observed for each 
carbon atom making it possible to follow deprotonation 
at each site in the molecule, In addition the range of 
chemical shifts is much larger than that observed in 
proton spectra making carbon-13 chemical shifts more 
sensitive to the state of protonation. Thus, it was 
thought that carbon-13 magnetic resonance spectroscopy 
might be a useful probe for the study of the acid-base 
chemistry of amino acids. The results obtained from a 
study of the acid-base chemistry of cysteine and related 
molecules by carbon-13 NMR are the subject of Part 2 of 


this thesis. 
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CHAPTER VII 


EXPERIMENTAL 


A. Chemicals 

D,L-=penicillamine, 2-aminoethanethiol hydro- 
chloride (Aldrich Chemical Company) and S-methylcysteine 
(Nutritional Biochemicals Corp.) were used as received. 
Cysteine (Nutritional Biochemicals Corp.) was recrys- 
tallized from water as the free base. 

All other chemicals were reagent grade and were 


used without further purification. 


B. Preparation of Solutions 

In general all solutions were prepared as des- 
‘eribed in Chapter II Section C, with the following 
modification to minimize ionic strength variations from 
sample to sample. Initially the pH was adjusted to ~13 
with potassium hydroxide or to (~lawith nitriciacid. 
Nitric acid or potassium hydroxide was then added and 
samples were withdrawn at the appropriate pH values. 
Using this procedure the ionic strength of a 0.2 M 
solution of cysteine for example varied from 0.9 at pH 


L360. 0.7 satu DHico. 


Ge pH Measurements 


pH measurements were made as described in 
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Chapter II Section D. 


D. Calculation of Hydrogen Ion Concentration 

In Chapter VIII the acid dissociation constants 
are determined as concentration constants. To convert 
pH meter readings, which are related to hydrogen ion 
activities, to PQH, the negative logarithm of the hydrogen 
ion concentration, the activity coefficient must be known. 
In order to calculate the activity coefficient for the 
hydrogen ion, the ionic strength must first be calculated 


using Equation (15). 


ie= VEIZ Cod. C15) 


where u represents the ionic strength of a solution con- 
taining i different ions of concentration C. and charge 
Zs In the solutions discussed in Chapter VIII the 
amount of ligand in the various charged forms was cal- 
culated assuming approximate values for the macroscopic 
constants. The concentrations of nitrate and potassium 
ions were calculated from the amounts of nitric acid and 
potassium hydroxide which were added to prepare the 
sample. The concentrations of hydrogen and hydroxyl ions 
were calculated from the solution pH. 

Once the ionic strength was calculated, the 


activity coefficient for the hydrogen ion was found using 
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the Davies Equation (18) 


log yo 0.511 (yn) 17? 


MIRE ey? TOR (yu) (16) 
n Tet as. Cu) ; 


where sas is-the activity coefficient of an ion in charge 
n. Using the meter reading as the negative logarithm of 
the activity of the hydrogen ion and the calculated 

activity coefficient, PH, the negative logarithm of the 


hydrogen ion concentration, was calculated. 


E. Non-Linear Least Squares Program 

In Chapter VIII the values for macroscopic 
dissociation constants were obtained from a non-linear 
least squares analysis of the chemical shift data. The 
non-linear least squares program used is based on the 
method described by Marquardt (19) and is available in 
the Chemistry Department of the University of Alberta 


as ENLLSQ. 


F. Carbon-13 Magnetic Resonance Measurements 


Carbon-13 magnetic resonance measurements were 
made as described in Chapter II Section F with the excep- 


tion that the number of pulses was reduced to 2K. 
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CHAPTER VIII 
CARBON-13 NMR STUDY OF THE ACID-BASE 


CHEMISTRY OF CYSTEINE AND RELATED MOLECULES 


x. Introduction 

In this chapter, the application of carbon-13 
NMR to the study of the acid-base chemistry of amino 
acids is investigated. The compounds chosen for this 
study are cysteine, which has been studied extensively 


by other techniques, penicillamine, 


- 


HSC (CH) ,CHCO, 


+ 
NH, 


and 2-aminoethanethiol 


HSCH.CH, 


op 


NH, 


First, the evaluation of macroscopic dissociation con- 
stants using carbon-13 chemical shift data will be des- 
cribed and then the attempts to evaluate microscopic dis- 


sociation constants will be discussed. 


B. Results 
1. Determination of Macroscopic Acid Dissociation 


Constants from Carbon-13 Chemical Shifts. 
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In this section, a method for the evaluation 
of overlapping macroscopic acid dissociation constants 
using carbon-13 chemical shift measurements is presented. 
The method is applied to the evaluation of macroscopic 
dissociation constants for the deprotonation of the amino 
and sulfhydryl groups of cysteine, penicillamine and 2- 
aminoethanethiol. The macroscopic dissociation constants 
for these deprotonation reactions will be described by K, 
and K, although in the case of Maminoethanethiol these 


3 


dissociations really represent Ky and K,- 


HoL 7 HL + H (17) 
rete fH 

Oy ers Seria (18) 

HL 2 H+ UL (19) 
2 EA 


In Equations (17)-(20) charges have been omitted 
for simplicity. 

Carbon-13 data for penicillamine, cysteine and 
2-aminoethanethiol were collected over the pH range 4-13 
where the amino and sulfhydryl deprotonations occur. 

This data is presented in Figures 30, 31 and 32. For 


cysteine and penicillamine, the assignment was done by 
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CHEMICAL SHIFT , ppm vs DIOXANE 
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Figure 30: pH dependence of the carbon-13 


chemical shifts of the five penicillamine 
carbon atoms in an aqueous solution con- 
taining 0.20 M penicillamine. 25°C. 
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Figure 31: pH dependence of the carbon-13 chemical 


shifts of the three cysteine carbon atoms in an 
aqueous solution containing 0.20 micvsteine. 25°C. 


: 
j 
i 
2 
Ka 
1 
— 4 
>. 
dine tA 
AOS ue nde 
t NE el pia 


er be 


CHEMICAL SHIFT ,ppm vs DIOXANE 


* 42.0 


46.0 


HSCHjCH NH» 


38.0 | 


PN Dec 8, SA Ra ST ke a 
400 6.00 8.00 10.00 2007 1260 


pes 


Figure 32: pHedependence of the carbon-13 chemical 
shifts of the two 2-aminoethanethiol carbon atoms 
in an aqueous solution containing 0.20 M 2-amino- 
ethanethiol hydrochloride. 25°C. oi 
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observing the off resonance decoupled carbon-13 spectrum, 
In 2-aminoethanethiol, the assignment was made by examining 
the magnitude of the carbon-13 shifts on deprotonation, 


Benesch and Benesch (6) have shown by ultraviolet measure- 


SH x NH,* SH 
ments that Ky Bey IN and thus that Ky = Ky - The 


assignment was made on the basis of the expected shift 
when NH,” deprotonates (20). One carbon resonance shifts 
2.16 ppm on the second deprotonation while the other 
shifts 5.24 ppm. The carbon resonance which shifts the 
most was assigned to the carbon 8 to the amino group. 

The observed shift for a particular carbon atom 
is the population weighted average of the shifts of that 


Carbon in forms: HaL, Hh and L;-that is 


2 

Sops = hay, Lae te Bae 4 Bay Se, 2 (21) 
where 
[HL] oe 

P Spe | ae ele = 
HL HoL ap HUBHIqlh se [hike 

re [HL] 
Pun i [HL] + [HL] + [L] (23) 

3 [1] 
Se ete sae [nae Da re 


The concentrations [HU], [HL] and [L] can be expressed 


in terms of the macroscopic acid dissociation constants: 
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eve to a try” (25) 
P, (HW + [H ]K, + KK, 
7 [H"]K, 
Pp = (26) 
HL +2 + 7 
{H } + [H K, + KK 
KK 
ers 
L +,2 + J 
ih.) + [H ]K, + KK 
and the observed shift can be expressed as 
+_2 + 
fe J °H DL + [H )K5 Spy, + KjK36,, 
5 Boy ped el a (28) 
OBS +_2 + 
[H ] + [H }K, + Kok 
Values for ou L and oy were obtained from the horizontal 
9 


portions of the shift versus pH curves at pH 4-6 and pH 
13-14 respectively. Values for [H’] and Sons were 
measured experimentally. Equation (28) was used to 


evaluate K K. and Sat, using a non-linear least squares 


2 3 


program. The values of K K, and Sur, obtained in this 


2! 
manner were then used to predict the observed shifts. 
The predicted and observed nites for cysteine, 2-amino- 
ethanethiol and penicillamine are given in Tables XIV, 
XV and XVI. Table XVII summarizes the results of these 
calculations and presents literature values of PK, and 
PK, for, comparison. | In most cases; ,it ist difficult: co 


make comparisons with the literature values for two 


reasons. First of all, the ionic strength used for our 
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experiments is considerably higher than that used in most 
of the experiments described in the literature. Secondly, 
the macroscopic constants reported here are'’in terms of 
concentrations whereas most of the literature values are 
mixed constants containing both activity and concentration 
terms. For purposes of comparison, a pH titration of 
cysteine was performed under conditions of ionic strength 
similar to that of the NMR experiment (uv = 0.8). Values 


Of 8.26 + O40) "andel0 se 2002 were obtained for pK. and 


2 
PK. These values compare well with those obtained by 


carbon-13 magnetic resonance. 


2% Determination of Microscopic Acid Dissociation Constants 
from Carbon-13 Chemical Shifts. 
At the molecular level, dissociation of the 
second and third protons from cysteine is described by the 


following scheme. 
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The observed shift is a population weighted 


average of the four forms, that is 


Sops ~ Proz * PyySzy t+ Prgrétrr * Pry ty (29) 
where 
[I] 
Py ll (3:0) 
[E] okpe ties [Ett] 2s. [Iv] 
[II] 
Prt = (3) 
[Ey Serre pIEn) 2 TLV] 
[III] 
ate Sth ALE ee ee riees (iv) 
[Iv] 
Pere (33) 
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The concentrations of species I - IV can be expressed 


in terms of the microscopic constants; 
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Keak 


h es ag We 
A hag Co udgyivds qudqutnemEciuecapse 
12 i3 Le igo Sy. 
and the shift can be expressed as 
: +,2 + | + 
H + 
suchas tig <5 © 5p ek i aii IRS Aas ee uni tas ee 
OBS + + + 
H + + 
He Bothy s TOTP) RD meek I SRT 45 
eee (38) 
which can be rearranged to the form 
Kan0 aed eS 6 
Sry ns 12 = : ITS2I EE fH] mt : : fH?) - 
ph Se IRS S32 
OBS 
k +k 1 
ane 12 | 13 tut} 4+ — [H']7 
a Tahal i i: Lo b32 
ecee (39) 


Equation (31) is of the form 


+ iy: 
5 baa SHSB MLO’ soho MO (40) 


Sel NES PMP a8 EE FS 
A non-linear least squares analysis of data obtained for 
a range of hydrogen ion concentrations will yield the 
coefficients a, b, c, d, and e. However, an infinite 
number of sets of Kio: Kj 3" Ki 32! Si7 and Osi will give 
the same b. In order to resolve the microscopic dis- 


sociation constants, values for Sir and Ssr7 are needed. 
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Equation (38) for the microscopic constants is 
similar to Equation (28) derived from the macroscopic 


constants. The shifts 6 and 6. are equivalent to 6 
HL L : £ 


and Sys The macroscopic product K K equals the micro- 


2 


scopic product k1 313° The coefficient of the [H"] term 


in Equation (28) is Koons while in Equation (38) it is 
ki o%rr as ki 3%rrr: Equation (38) can be shown to be 
exactly equivalent to Equation (28) by the above replace- 


ments and the following: 


San = Prr8sr * PrrrSzrr: (41) 
ee nes eRe a oe 
SSPE ER [Vea ae o 
HL Ky K, 
KoSnr = ¥yoStr * k13%rrr- (43) 


An alternative way of obtaining values for the 
microscopic constants is to use the value Sur, which is 


obtained for the shift of the monoprotonated species 


from the evaluation of the macroscopic parameters. Since 
oH Pte are 4 TT ETT be?) 
where Fry and For represent the fractional concentrations 


of species II and III and 


F ce al, (45) 
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be 
Pee meta 3! = ep (46) 
HEAL ete 
Also, 
Tie eae 
Fa ee cats eae 
ALALAL 3 
Therefore 
a: 2 
OO ee aT aa ew oe) 
AGALAL 2 LZ 
and 
pee (49) 


Thus jaar t Orr and Orr are known, the microscopic constants 
can be evaluated directly from Our: Attempts were made 
to determine Oot and Orrr using model compounds. 
S-methylcysteine was the first model compound 
used... At. pH 15.32, the, shifts: ofsthe Cyr Ce and “cH, carbons 
ALC. SUG uO wp, FO 2:04. a7, ana 51.95 ppm respectively, while at 
pHa 2. Av sthescorrespongding shi ftisiware sl 655,02 735 . and 
Sle opp welhus. the Cy resonance. shistts wy ede. Sirs ppm 


while the C, resonance shifted by 4.69 ppm. If one 


B 
attributes these shifts to the deprotonation of the amino 
group and assumes that the Cy and Ce resonances of cysteine 


shift by the same amount when the amino group is de- 
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protonated, the shifts of species II and III may be pre- 
dicted. For example, the Cy resonance of cysteine is at 
10.59 ppm in species I and 6.20 ppm in species IV, giving 
a total shift of 4.39 ppm. If 1.51 ppm of this shift 

is due to the amino deprotonation, the remainder, 2.88 
ppm, must be due to the sulfhydryl deprotonation. Thus 
the oF Shai tein ispecies Meise) 059-288. 70r. 7.71, ppm: 
In species III, the Cy Shafts bO.50 <1 25) | or 79,08 
ppm. The shift Son for the C, carbon was found from the 
non-linear) least squares: fit to bel8./72 .ppm,.. One.obtains 
a value of 0.263 ‘for Fir from Equation (46) and the above 
shifts, 


A similar analysis of the C, chemical shifts gave 


B 


a valie of (0,920 for For This value doesn't agree with 


the value obtained from the Cy data nor with the known 


relative acidities of the sulfhydryl and amino groups. 


Attempts were jmade- to perform similar calculations 


using glutathione as the model compound. The cysteinyl 
residue of glutathione can be considered an N-substituted 
cysteine and the shifts of ene CYS-C, and CYS-Cp carbon 
resonances attributed to sulfhydryl deprotonation. 
The shifts of the intermediate species, II and III, can 
be calculated using the observed glutathione shifts by 
the procedure described above for S-methylcysteine.., By 
BLE 


treating the Ce data, in’ this manner, F was found to be 


0.137, while a value of 0.748 was obtained from the 
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analysis of the Ce nea Again the agreement is unsatis- 
Factory: 

Benesch and Benesch (6) concluded from ultra- 
violet measurements that the first deprotonation in 2- 
aminoethanethiol is from the sulfhydryl group and the 
second from the amino group. They found that there was no 
Overlap or mixing of the deprotonation reactions. MThere- 
fore 2-aminoethanethiol can be considered as a model 
compound and its shifts used to calculate Say and Ssrr 
for cysteine. Since the total shift on deprotonation of 
both groups is not the same in 2-aminoethanethiol as it 
is in cysteine, the shifts due to deprotonation of a 
specific group were estimated by assuming that the fraction 
of the total shift due to deprotonation of the sulfhydryl 
group is the same in 2-aminoethanethiol as it is in cysteine. 
Analyzing the data in this way, one obtains a value of 
Or29 Lor fF. from the C, results. The Cy data showed 


TL 8B 
Orr and Ss3r both to be less than Our, so that it was 
impossible to obtain Fir from the C. shatt data, 
The results of the above calculations of Sar 


and Ort for cysteine from the shifts observed in model 


I 
compounds are summarized in Table XVIII. The values 


predicted for the ratio R of the sulfhydryl-deprotonated 
form to the amino-deprotonated form are also included in 
this table. To estimate the approximate value of Re pH 


titrations of cysteine and S-methylcysteine were performed 
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Table XVIII, Predicted Values of Sit and Ssar using 


Model Compounds. 


Model Compound Carbon 
S-methylcysteine Ch 

Cg 
Glutathione ae 

Cg 


2-aminoethanethiol 2: 


a : 
ppm versus dioxane 


PRED 


oT 


CLR Bal 


39,248 


7.945 


40.045 


8.545 


40.205 
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: ratio of sulfhydryl-deprotonated intermediate to amino- 


deprotonated intermediate. 
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at an ionic strength of 0,8 which is approximately the 
ionic strength used in the NMR experiments, Macroscopic 
constants for PK, and PK, of cysteine were determined to 
be 8.26 and 10.38. For S-methylcysteine a value for pK, 
of 8.80 was obtained. By assigning this value to the 
microscopic constant Pk) 3 of cysteine, the microscopic 
constant Pk, > was determined to be 8.41. The value of R 
obtained from these data is 2.48. 


C. Discussion 
1. Evaluation of Macroscopic Constants From Carbon-13 
Chemical Shift Data. 

The macroscopic constants obtained for cysteine, 
penicillamine ane 2-aminoethanethiol from carbon-13 
chemical shifts indicate that this is a reliable method 
for determining macroscopic dissociation constants even 
when the deprotonation reactions overlap. The results 
listed in Table XVII compare favorably with those avail- 


able in the literature. 


2. Evaluation of Microscopic Dissociation Constants From 
Carbon-13 Chemical Shift Data. 
The evaluation methods described earlier in 
this chapter indicate that tnere*arerdittscultlres™an “using 
carbon-13 chemical shift data as a means of obtaining 


microscopic dissociation constants when the chemical shift 
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of the carbon atom is affected by the deprotonation of 
both groups. In such cases, it is necessary to know the 
shifts of the two monoprotonated isomers in order to 
determine the microscopic constants. These shifts cannot 
be obtained directly since both isomers are always present. 

One approach which was tried was the use of 
model compounds to estimate the shifts of the mono- 
protonated forms. This approach was not successful, which 
is not too surprising in view of the assumptions involved 
in the model compound approach. First of ali, carbon-13 
shifts are known to be extremely sensitive to confor- 
mational changes (26). It is not possible to take into 
account the differences in the conformations of the model 
and the actual compound nor the conformational changes 
that occur on deprotonation, 

Another assumption involved in the model compound 
approach is that the shift of a given carbon resonance 
upon deprotonation of a certain functional group is always 
the same no matter what the state of other functional 
groups. In other words the meted Aarne Lo that, tora, par- 
ticular carbon atom, Sor - oy equals Siy - Sort 

Recently Walters and Leyden (15) used the model 
compound approach in a determination of the microscopic 
ionization scheme of cysteine using pmr spectroscopy. 

Not surprisingly, the results they obtained are in dis- 


agreement with the results obtained by other workers using 
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more direct methods (6,7). Both CMR and PMR methods have 
the same basic problem: the determination of the shifts 
of the m ec i : 
e€ monoprotonated species, Sir and Sirtr 
Sion which can be derived to describe the observed shift 


Any expres- 


as a function of the microscopic constants necessarily 
involves oor and Osrtr Until good estimates of these 
parameters can be obtained, it is impossible to evaluate 


the microscopic ionization constants from NMR data. 


3, Conformational Dependence of Carbon-13 Shifts. 

In the results for penicillamine presented 
in Figure 30, one of the methyl carbons shifts first in 
One direction as one proton is removed and then in the 
other direction as the second pEeeon is removed. Pre- 
sumably this occurs because the major conformational form 
is different in the deprotonated species from that in 
the protonated species. 

TLan-Dinh et.al (20) “have recently applied 
carbon-13 magnetic resonance to the problem of determining 
Con romans of the four amino acids: alanine, valine, 
leucine and isoleucine. Their results for valine are 
representative and are discussed below. The three con- 


formers of valine can be written as 


CH CH, 


- + ri + 
= NH,* 0,¢ NH 2G NH 
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The position of the carbon-13 resonance for each 
methyl carbon changes uniquely with pH. One methyl 
resonance shifts as both the carboxyl and amino groups are 
deprotonated, while the other methyl resonance remains 
constant. From this it was concluded that one methyl carbon 
must be localized in the <coo , NH, "> ZOne aS in conformer 
I or II and that the relative proportions of conformers 


I and II do not vary with pH, 


In penicillamine the three conformers are: 


: ‘iby 
ee SH She CH, CH; he SH 
once : has oe ae 
ah BAe CH, 


The anomalous shifts shown by the two methyl groups (Figure 
30) can be explained in terms of conformational changes 
as deprotonation proceeds. On going from H,L to HL both 
methyl groups shift by a similar amount which indicates 

+ 


that both spend some time in the critical <Ccoo , NH, > 
zone. Thus a mixture of conformers II and III exists. 

On going from HL to eee a change in the conformational 
population occurs. Only one methyl resonance shifts 
further downfield on deprotonation suggesting one con- 
formation is predominating, most likely III since in this 


rotamer the two negatively charged groups, COO and S_ 


are trans to each other. 
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